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Abstract The ground states of all even-even nuclei have angular momentum, /, 
equal to zero, / = 0, and positive parity, vr = +. This feature was believed to be a 
consequence of the attractive short-range interaction between nucleons. However, in 
the presence of two-body random interactions, the predominance of I'^ = 0"*" ground 
states (0 g.s.) was found to be robust both for bosons and for an even number of 
fermions. For simple systems, such as d bosons, sp bosons, sd bosons, and a few 
fermions in single-j shells for small j, there are a few approaches to predict and/or 
explain spin I ground state (/ g.s.) probabilities. An empirical approach to predict 
/ g.s. probabilities is available for general cases, such as fermions in a single-j 
{j > 7/2) or many-_7 shells and various boson systems, but a more fundamental 
understanding of the robustness of g.s. dominance is still out of reach. Further 
interesting results are also reviewed concerning other robust phenomena of many- 
body systems in the presence of random two-body interactions, such as the odd-even 
staggering of binding energies, generic collectivity, the behavior of average energies, 
correlations, and regularities of many-body systems interacting by a displaced two- 
body random ensemble. 

PACS: 05.30.Fk, 05.45.-a, 21.60Cs, 24.60.Lz 

key words: I g.s. probabilities, g.s. dominance, random interactions, correla- 
tion, collective motion, average energies. 
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1 Introduction 



Atomic nuclei provide ideal laboratories to study the features of microscopic many- 
body systems with finite number of constituents less than ~ 300. They are complex 
objects with many degrees of freedom and exhibit almost all features found in other 
many-body systems. The regularities in atomic nuclei in the presence of random two- 
body interactions therefore provide an excellent window to study general features 
of low-lying states of many-body systems. Therefore, although discussions in this 
article are originated from nuclear structure studies, the results can have many 
implications for other fields as well. 

The Gaussian orthogonal ensemble of random matrices was first proposed by 
Eugene Wigner in Ref . [1] , which was a revolutionary thought in understanding the 
spacings of levels observed in resonances in slow-neutron scattering on heavy nuclei. 
The two-body random ensemble (TBRE), which will be used in most examples in 
this article, was introduced to study statistical properties of spectra of many-body 
systems, by French and Wong in Rcf. [2], and by Bohigas and Florcs in Ref. [3]. 
Rcf. [4] presented a self-contained account of random matrix physics in quantum 
systems concerning spectrum and strength fluctuations. The latest review along 
similar lines is presented in Ref. [5]. 

We also note that there recently appeared other review articles on random ma- 
trix ensembles for finite particle systems but with different focuses. In Ref. [6], the 
focus was on the links between many-body pairing, as it evolves from the original 
nucleon-nucleon force, the manifestations of superfiuity in nuclear matter (say, neu- 
tron stars), and pairing in atomic nuclei. In Refs. [7, 8] the discussion was focused 
on statistical mechanics and onsets of chaos in finite many-body systems. In Ref. 
[9], the focus was on the geometric chaoticity of angular momentum couplings, its 
possible implications on the energy ccntroids, and the multipole collectivity in the 
presence of random two-body interactions. In this article we shall focus on the fea- 
tures (particularly, orders and correlations) of low-lying levels of many-body systems 
in the presence of random two-body interactions, explaining both the observations 
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and the present status towards "understanding" these features. 

In many-body systems such as molecules and atomic nuclei, the interactions by 
themselves have no trace of symmetry groups for vibrational or rotational modes. 
However, the low-lying states often exhibit a pattern suggestive of symmetries for 
these modes. One may ask to what extent the low-lying states acquire order from 
the basic properties of interactions such as rotational invariance and possibly other 
symmetries such as isospin invariance. In other words, some properties such as 
vibration or rotation might dominantly occur in the low-lying states of many-body 
systems while the others might occur only with small probabilities. 

Atomic nuclei with an even proton number Z and an even neutron number N are 
examples which follow these lines. The angular momenta and parities {I'^) of the 
ground states of even-even nuclei are always 0^, and the structure of the low- lying 
states is characterized by a tripartite classification [10, 11], i.e., seniority region, 
anharmonic vibrational region and rotational region. Also there exists an odd-even 
staggering of binding energies, etc. One can ask whether these features are robust 
for general many-body systems. This can be studied by permitting interactions to 
be more and more arbitrary. 

This question of robustness was first studied by Johnson, Bertsch and Dean in 
Ref. [12], where a dominance of — 0+ ground states (0 g.s.) was obtained by 
using the TBRE Hamiltonian. In Ref. [13], the g.s. dominance was found to be 
insensitive to the monopole pairing interaction, and to be related to a reminiscence 
of generalized seniority suggested in Refs. [14, 15]. In Refs. [16, 17], it was found 
that the sd bosons produce both vibrational and rotational spectra as well as the 
g.s. dominance in presence of the TBRE Hamiltonian. In Ref. [18], it was found that 
odd-even staggering of binding energies in finite metallic grains and metal clusters 
[19] arises from purely random two-body interactions. These interesting results 
suggest that the above features are not only the consequences of attractive pairing 
interactions, but are much more general than they were previously assumed. Many 
authors made efforts to seek the origin of these observations and to look for more 
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robust features along these lines. This article aims at reviewing these achievements. 

As was emphasized by Fcshbach [20], the studies of atomic nuclei provide us with 
many "universals" of the physical world. In this article we will show that microscopic 
many-body systems in the presence of random interactions provide us with a new 
method to discover and to study universal features of microscopic systems in nature. 
Discoveries and understandings of these interesting patterns are very exciting topics 
in physics. 

In Sec. 2 we shall define the Hamiltonians for systems such as nucleons and/or 
other fermions in a single- j or many-j shells, bosons with single spin I or many Vs. 
In the process of doing this, we shall also give a brief introduction to the nuclear 
models at a basic level, which will be helpful to those who are not famihar with the 
nuclear structure theory. We then define the two-body random ensemble that we 
take in this paper. 

In Sec. 3 we shall concentrate on statistics of the distribution of spin / in the 
ground states in the presence of the TBRE Hamiltonian. We shall not restrict 
the discussions to the g.s. probability (denoted as -P(O)) which is obtained by 
diagonalizing the TBRE Hamiltonian, but also study other I g.s. probabihties 
(denoted as P{I))- We shall go to systems with odd numbers of fermions to study 
other I g.s. probabilities and parity distributions of the ground states as well. 

In Sec. 4 we shall review the efforts to understand the g.s. dominance of simple 
systems suggested in earlier works, where one can evaluate / g.s. probabilities. The 
systems are restricted to d-, sp- and sci-boson systems, and fermions in a single-j 
shell with j < 7/2. 

In Sec. 5 we shall go to more complicated systems, such as fermions in a singlc-j 
(for large j) or many-j shells and sdg bosons, etc., for which the understanding of 
P(/)'s is not yet available. However, an empirical approach will be suggested to 
predict -P(/)'s of these systems. This empirical approach also shows that the g.s. 
dominance is related to certain interactions with specific features. Some properties 
of /max g-S. probabilities will be found and explained, and some features of the energy 



7 



gaps for g.s. and the n-body matrix elements of / = states will be pointed out. 

In Sec. 6 wc shall discuss other features of many-body systems interacting by 
random interactions, including average energies, collectivity, and normal ordering of 
spin / for the yrast states. We shall show that the energy centroids of spin / states 
(denoted by Ej) have large probability to be the lowest energy only if / ~ /min 
or 7 ~ /max- Compact trajectories of £'/'s can be found. Bijker and Frank [16, 17] 
discovered that both the vibration and the rotation are robust in the low-lying states 
of sp- and sd-boson space. For fermionic degrees of freedom, say, nucleons in the sd 
shell, rotational spectra do not appear if one uses a general two-body Hamiltonian. 
We will show that the quadrupole-quadrupole correlation is very important to get 
a collective rotational spectrum. We will also see in Sec. 6 that spin I in the yrast 
bands obtained by the TBRE Hamiltonian has a large possibility to be normally 
ordered (i.e., the sequence 0, 2, 4 ■ ■ ■). In the fourth part of Sec. 6 we shall review 
results of constraints of the TBRE Hamiltonian in order to reproduce some global 
features of atomic nuclei. 

In Sec. 7 we summarize this paper and conclude on this interdisciplinary field. 

In the Appendix we listed mathematical notations of this paper. 
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2 Space and ensembles of the Hamiltonian 



In this Section we define the Hamiltonians to estabhsh notations and conventions 
used in this paper. Towards this goal, we present a brief introduction to the nuclear 
shell model [21, 22, 23] and its boson approximation, the interacting boson model 
(IBM) [29, 30, 31, 32]. This helps those who are not experts at the nuclear structure 
theory, which is frequently involved in the studies of the present subject, to catch 
easily the essential point of physics. More detailed discussions can be found in Ref. 
[15]. We also explain the random samplings and the two-body random ensemble. 
These are very basic concepts to study, e.g., the question how large a probability is 
for a certain spin I to be the spin of the ground state. 

2.1 Hamiltonian of the shell model and its approximations 

The simple picture that particles move independently in a one-body potential well, 
which well describes the atoms, metals etc., is also a key to the nuclear structure 
theory. Here the one-body potential is given by the mean field which is produced 
by all nucleons except the "single" nucleon in consideration. The Hartree-Fock 
method provides us with an approach to derive a single-particle potential from the 
two-body interactions, with the requirement that the energy for A-particle Slater 
determinants should be a minimum. Numerical results of Hartree-Fock calculations 
indicate that the harmonic-oscillator wave functions are good approximations to 
the self-consistent wave functions. A breakthrough was made by Mayer and Jensen 
[21, 22, 23] through introducing the strong spin-orbit coupling in the nuclear force. 
This nuclear shell model (SM) makes understandable a vast amount of experimental 
data of spins, magnetic moments, isomeric states, and the "magic numbers" for 
proton number Z or neutron number N equals 2, 8, 20, 28, 50, 82, 126. The SM 
has been playing a central role in understanding the nuclear structure since it was 
established in 1949. It is defined by a set of spin-orbit coupled single-particle states 
with quantum numbers rirljm, corresponding to the radial quantum number (n^), 
orbital angular momentum (l), total angular momentum (j) and its 2;-component 
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(m), respectively. The parity for each orbit is determined by (— )^ Because the 
Hamiltonian is rotationally invariant, the single particle energies are independent 
of m. For all single-particle states within one major shell, the value of each j 
is unique, so one usually suppresses the quantum numbers nl and uses only the 
total angular momentum j to label a certain single-particle level. When we define 
the Hamiltonian below, we shall use Ji,j2,j3,j4 to label different single-particle 
levels. In atomic nuclei there are two types of particles, neutrons and protons. 
One introduces the isospin quantum number t, which is mathematically very similar 
to spin s, with t — 1/2 and rrit — ±1/2 corresponding to protons and neutrons, 
respectively. Because t for both protons and neutrons is the same, one can keep 
only rrit to label the single-particles states. The single-particle energies for protons 
and neutrons are different, with the proton ones higher mainly due to the Coulomb 
force. It is usually believed that at least one major shell is necessary to adequately 
describe low-lying states of a given nucleus. 

One property of nuclear systems is that the residual interaction between the 
valence nucleons, i.e., the part of the interaction which can not be covered by the 
mean-field part, is very strong and leads to strong configuration mixings. Therefore, 
one has to diagonahze the Hamiltonian in the shell model space to obtain the "true" 
eigenvalues and wave functions. 

The shell model Hamiltonian that is usually taken includes a one-body term 



-^1 ~ ^ ^jmt^\m,mt^jm,mtj (l) 



jmmt 



and a two-body term 



^2 = i V ( 1 + '^iU2 ) ( 1 + ) GjT (jl J2 , J3J4) 

jlj2j3j4,JT 

E ^Uiii2)a,^03j4)a, , (2) 



where 



MjMt 



'2 2 ^ ' / Jimi, mtj "'J2m2,mt2 ' 
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In Eq. (3), (jimi,j2m2|JMj) denotes the Clebsch-Gordan coefficient. In Eqs. (2-3) 
J (T) is the total spin (isospin) resulting from one particle in the j'l orbit and another 
one in the j2 orbit, or one particle in the js orbit and another one in the orbit. 
Mj and Mj^ are the z components of J and T, respectively. The a^^ and 
single-particle creation operators. GjT{jij2, jsj^ys are two-body matrix elements 
defined by 

Gjt(JiJ2, J3J4) = {{jlj2)MMT 1^1 {j^ji)MMT)^ (4) 

for a two-body interaction, V. The GjT{jij2, jsjiYs are independent of M and Mr- 
The Hamiltonian of Eq. (2) respects parity conservation, which means that the 
parity product for the four orbits J1J2J3J4 is positive. The Hamiltonian Eq. (2) also 
respects rotational and time-reversal invariance, which means that V in Eq. (4) is 
a scalar and takes real values. 

In this paper we shall also consider systems in which there are only one type 
of valence particles (i.e., neutrons or protons). For a single-j shell one then can 
suppress the labels jij2j3j4 (ji = J2 = ia = = j) and T {T = 1) and simply 
denote the two-body matrix elements by Gj, where the allowed values of J are 
J—0, 2, • • •, 2j — 1. The number of G/s of fermions in a single-j shell is j + |. 
It is noted that the "one type of particle systems" here should be understood in a 
broader sense: the identical particles are not necessarily identical valence nucleons, 
they can be also other types of identical particles (say, electrons). Therefore, we 
refer to these systems as "fermions" in a single-j shell or many-j shells. 

In this paper we shall also use a restricted separable Hamiltonian for fermions 
in many-j shells: 

H^H, + Hp + Vp^, (5) 

where Hi, Hp and Vph arc the spherical single-particle energy term, generalized 
pairing, and particle-hole type interactions, respectively. The definition of Hp is as 
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follows, 

Hp = Vo + V2 + --- . (6) 

Here 

i 

M jlj2 

with 

^<°> = (A*'"')' , ^S' = {AW'f , ,uuy = . (8) 

The particle-hole interaction Vph takes the form 

Vpi,^ kY^QmQm ^ , (9) 

M 

where 

Qm= E q{jij2){-)"{jirn,,j2rn2\2M)al^^aj,^rn, ■ (10) 

jij2mim2 

The operator Q is called the quadrupolc operator, and the interaction J2m QmQm is 
called the quadrupole-quadrupole interaction. The validity of truncating two-body 
interactions to the pairing and the quadrupole-quadrupole interaction in nuclear 
physics was studied by Kumar and Baranger [24, 25]. 

The difficulty to apply the nuclear shell model is that the dimension of the con- 
figuration space involved in a major shell increases very rapidly with the number 
of valence nucleons. For medium and heavy nuclei, the number of states is usually 
too huge to handle even for a very good computer. The status of diagonahzation of 
the shell model Hamiltonian was reviewed in Ref . [26] . Because of this difficulty, 
many efforts were made to simplify the shell model approach. Important milestones 
include the Bohr-Mottelson-Rainwatcr collective model (or the so-called geometric 
model) [27, 28] and the interacting boson model (IBM) [29, 30, 31, 32]. The geo- 
metric model was introduced by Bohr, Mottelson, and Rainwater in 1950-1953 as an 
alternative to the shell model. The link between the geometric model and the shell 
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model was established by Elliott [33] in 1958. Intrinsic deformation is introduced in 
the geometric collective model. In the IBM the building blocks of the Hamiltonian 
and the model space are s bosons with spin zero and d bosons with spin two, which 
are interpreted as correlated S nucleon pairs with spin zero and D nucleon pairs 
with spin two. Namely, the IBM is a phenomenological model which first truncates 
the full shell model space to the S and D pair subspace and next maps the S and D 
pair subspace to the s and d boson subspace. The total boson number is conserved 
in the Hamiltonian and equal to half of the valence nucleon number for an even-even 
nucleus. The parity of s and d bosons is positive. The relationships between the 
IBM and the geometric description were discussed in Refs. [32, 34, 35, 36]. 

The IBM in which the distinction between protons and neutrons is left out is 
called the IBM-1. The IBM-1 Hamiltonian [32] that we use in this paper is as 
follows, 

m L=" ~ 

1 v-r/,+ ,+\(2) 



2 

L=0,2,4 ^ M 



IM 



Y.\(d^d^yUMS + h.c. 



+ ^ssss 

S^S^SS , (11) 



where 



and we take this notation for bosons hereafter. 

Similarly, we shall also go to the sdg and sp systems. The spinP^'*^ for a g boson 
and a p boson is 4+ and 1^, respectively. The sdg boson systems have been studied 
by including g bosons as well as sd bosons in describing deformed nuclei. The sdg 
Hamiltonian is very similar to that of the sd lBM-1 Hamiltonian except that there 
is another one-body term for the g bosons and there are many more two-body terms 
(totally 32 two-body terms) in which g bosons are involved. One sees Ref. [32] for 
details of the sdg IBM. The sp interacting boson model is also called the "vibron" 
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model [37], which is used to describe both rotations and vibrations of the diatomic 
molecules. The total parity for a sp-boson system is given by {—)^, where / is 
the total spin of a certain state for this sp-boson system. Like the sd IBM, the 
total number of bosons is also conserved in the sp IBM (or the vibron model). The 
Hamiltonian for the sp bosons is as follows, 

1 f / . a(o) 



L=0,2 M ^ M 



+ 2 ^PPss 



We shall also use a schematic sp-boson Hamiltonian 

H = -cosx J^PmP"^ + 4/^^-1 N (5^5^ - pV) {ss - pp)^^^ , (13) 

which contains the basic features of the model, to exemplify the so-called mean-field 
approach to spin I ground state probabilities. 

The success of the sd IBM in describing the low-lying excitations of atomic 
nuclei stimulated calculations within the subspace constructed by correlated S and 
D pairs. Similar to the IBM, Ginocchio constructed symmetry-dictated SD pairs 
[38]. This approach was further studied by Wu et al. and called the Fermion 
dynamical symmetry model [39]. To overcome the restriction in Ginocchio's model 
that the structure of the SD pairs is very specific, Chen [40] developed recursion 
formulas which are apphcable for arbitrary pairs. Based on this technique, a nucleon 
pair approximation (NPA) of the shell model was suggested in Refs. [41, 42). The 
NPA Hamiltonian in this paper is defined by Eq. (5). 

Parallel to fermions in a single- j shell, we shall also study boson systems with 
spin I. The Hamiltonian for a boson system with spin I is written as 

evenL 



where h\ is the creation operator for bosons with spin /. The number of two-body 
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matrix elements for a boson system with spin / is / -|- 1, and {hibi)^M^ = i{^\b\^ *^ 



In the numerical calculations of this paper, the one-body terms in the above 
Hamiltonian will be neglected. It is noted that the inclusion of one-body terms does 
not change the statistics substantially. 

2.2 Monte Carlo samplings 

The name "Monte Carlo" arises from the "random" character of the method and the 
famous casino in Monaco. The Monte Carlo method is well known and a powerful 
method in almost all fields of physics. For instance, there are 10 review papers with 
"Monte Carlo" in their titles in the journal Physics Reports. This method has been 
applied to many problems, such as the evaluation of high-dimensional integrals, the 
shell model Monte Carlo method in nuclear physics [26], the quantum Monte Carlo 
simulation in solid state physics [43], etc. 

It is usually difficult to investigate questions such as how large the probability 
is for a system to have a rotational spectrum or to have a certain spin / ground 
state for the full space of two-body interactions. In particular, when the number of 
two-body interactions is very large, it is impossible to evaluate those probabilities 
by summing over the results obtained from the usual quadrature points throughout 
the full parameter space. 

The approach by random samplings of the ensemble is a very efficient and simple 
way to study the regularities of low-lying states of many-body systems throughout 
the possible parameterizations. The advantage by Monte Carlo samplings can be 
seen, for example, in the calculation of integrals with many variables. Suppose one 
calculates an integral with D variables, then the accuracy by the usual quadrature 
method is p~^/^ while that of the Monte Carlo calculation is Here p is the 

number of quadrature points of the conventional quadrature method, or the number 
of samplings of the Monte Carlo method (obviously, p » 1). When the integral is 
of very high dimension (i.e., D >> 1), the Monte Carlo samphng is an efficient way 
to evaluate the integral, while the usual quadrature method becomes completely 
impossible. 
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There are many computer programs to generate random numbers with specific 
distributions. Of course, a sequence of numbers generated by a deterministic pro- 
gram can never be truly random. Nevertheless, the sequence looks "random" (called 
pseudorandom) enough for practical purposes, and thus can be used as if they were 
truly random. A good discussion of Tiniform random number generators and of tests 
whether or not these "random" numbers work properly can be found in Ref . [44] . 

In studying the regularities of many-body systems in the presence of random 
interactions, one thousand samplings present very stable statistics and thus are 
deemed to be "reliable" samphngs. One then says that a "rehable" statistics is 
obtained for aspects such as spin I ground state probabilities. 

2.3 Two-body random ensemble 

The two-body random ensemble is defined as follows. The two-body matrix elements 
are independent of each other, and take random values which follow Gaussian dis- 
tribution with an average being zero. The distribution width is set to be 1 for 
diagonal two-body matrix elements and 1 / y/2 for off-diagonal two-body matrix el- 
ements. More specifically, if Gjr0'ij2, ^3^4) and Gj>T'{j'ij2,j'2,j'i) are two arbitrary 
(different) two-body matrix elements for valence nucleons in an open shell, one has 

/. . . . XX 1 / [G'jT(jlj2,j3i4)]\ 
P[GjT{jlj2,j3j4)) = ^==CXp( — ), 

^(r< (A' A' A' A'W ^ ^^^^( \GjiT' (i^i2^i3i4)]^ ^ /i r\ 

V 2'nx' 2,x' 

where 

^ ^ r 1 if \{h32)jT) = \{hk)jT) ^ r 1 if \{mj'T') = uf.yr) 

y \ otherwise ' 1 | otherwise 

One sees for the ensemble average that 

i2\ 



(16) 



{GjT{3i32,3zk) Gj'T'{j[j2,333'd) = ■ 
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A two-body random ensemble such defined is called the "TBRE". In this paper 
wc shall discuss fermions in a singlc-j shell or many-j shells, and bosons with spin I 
or more than one spin. The Hamiltonians for these systems are defined in subsection 
2.1. The TBRE Hamiltonian for each of these systems is given by an ensemble in 
which two-body matrix elements, such as Gj's for fermions in a single-j shell or G^s 
in Eq. (14) for bosons with spin I, or e's and c^'s in Eqs. (11)-(12), are taken to be 
random numbers following the distribution of Eq. (15). The one-body interaction 
parameters are set to be zero. 

In Ref. [12], Johnson et al. also used an ensemble for which the width x has 
a J- and T-dependence: the width equals to the above x multiplied by a factor of 
(2j+iK2r+i) • This was called the random quasiparticle ensemble (RQE). In Refs. 
[9, 45, 46] the MSU group took two-body matrix elements to be random numbers 
which are uniformly distributed between —1 and 1. 

It is noted that the properties obtained by using the TBRE, the RQE and the 
uniformly distributed random two-body matrix elements are essentially similar to 
each other, although the statistics might be somewhat different. This was already 
noted in the early paper by French et al. [2] for the spectral statistics. 

The results of this paper are based on 1000 samplings of the TBRE Hamiltonian. 

2.4 Summary of this Section 

In this Section we have defined the Hamiltonian corresponding to valence nucleons 
in many-j shells, fermions in a single-j shell or many-j shells, sp bosons, sd bosons, 
sdg bosons, and bosons with spin I. We have presented the connection between the 
Hamiltonians of these different systems by reviewing the nuclear shell model and its 
various approximations: the Bohr-Mottelson geometric model, the Arima-Iachello 
interacting boson model, etc. This knowledge is helpful in estabhshing notations 
and conventions in the paper. The separable Hamiltonian is taken for the SD pair 
approximation of the shell model. The sp boson model was introduced to describe 
the rotation and vibration of diatomic molecules, but the spirit is very similar to 
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the interacting boson model. 

We have introduced very briefly the concept of Monte Carlo samphngs. The 
Monte Carlo approach is widely used in various subjects. The term "Monte Carlo" 
appeared in titles of hundreds of monographs, and Monte Carlo techniques were 
used in numerous papers. Here we have explained why the Monte Carlo samplings 
are very efficient in studying the regularities of many-body systems in the presence 
of random interactions. We point out that 1000 sets of random interactions produce 
statistics which is "accurate" enough for our discussions. 

We have also defined the "two-body random ensemble" (TBRE). One-body in- 
teraction parameters in the Hamiltonian are set to be zero, and two-body interaction 
parameters are independent of each other and chosen to follow a Gaussian distribu- 
tion of width for the off-diagonal terms and 1 for the diagonal terms. We note 
that some authors took non-Gaussian distributions, but the pattern exhibited by 
the statistics such obtained is very similar. 
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3 Regularities of the ground states 



In this Section we shall go first to two simple cases: fermions in a single-j shell and 
bosons with spin I. We take four particle systems for these two cases and change 
j and I values, and focus on spin I ground state probabihties based on the TBRE 
Hamiltonian. Then we review the results for more complicated systems based on 
both the TBRE Hamiltonian and displaced TBRE Hamiltonian. We also study 
the regularity of parity distribution in the ground states for nuclei with random 
interactions and the effect from higher rank Hamiltonians. 

3.1 Fermions in a single- j shell 

Before going to complicated cases, let us first come to a few simple examples. A 
system with four fermions in a single-j shell is the simplest but non-trivial case 
for fermionic degrees of freedom. This case was treated by using both the TBRE 
Hamiltonian in Ref. [47], and a uniformly distributed random two-body Hamiltonian 
assuming values between —1 and 1 in Ref. [45]. The results obtained by using these 
two ensembles are essentially the same although the distributions of the ensembles 
are very different. Fig. 1 shows the results of a few important / g.s. probabilities 
calculated by the TBRE Hamiltonian, with j ranging from j — 7/2 to j — 33/2. 
Some of the principal outcomes are: 

1) The probabihty of ground states with odd / is much smaller than that of their 
neighboring even values of spin, even if the corresponding numbers of states in the 

configuration are comparably large. 

2) The unique /max state has a large probability to be the ground state, although 
this probability decreases with j ^. 

3) The 7 = 2"'" and 4"*" states have large probabilities to be the ground state. This 

indicates that small and even angular momentum states are favored as the ground 

^This holds also for the case of random intt^actions which distribute uniformly between —1 and 
1. The results of Ref. [45] which claim that this probability staggers rapidly and becomes for 
several single j are erroneous. 
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states for an even number of fermions in single-j shells. 

4) The g.s. dominance obtained by using the TBRE Hamiltonian is not a 
"rule" without exceptions. In Fig. 1, P(0)'s of four nucleons in a j = | and j = y 
shell are smaller than the corresponding P(2)'s. 

5) An interesting oscillation of P(0)'s with respect to j was noticed in Refs. 
[46, 47]. This oscillation is synchronous to an increase of the number of / = 
states. Namely, P(0) staggers when {j — |) = 3/c {k is a positive integer), coinciding 
with an increase of the number of 7=0 states which is given by the largest integer 
not exceeding (j - |)/3 [48, 49]. 

Because the g.s. is dominant for four fermions in a single-j shell with the only 
two exceptions j =7/2 and 13/2, one easily gets the intuition that the P(j)'s might 
be large for five fermions in a single-j shell. This was found to be indeed the case 
in Refs. [45, 47], although the P(j)'s are not as large as the P(0)'s. Fig. 2 presents 
a few examples for n = 5 in which P(j)'s are indeed very large in comparison with 
other P(/)'s. 

One connection between the dominant P(0) for n — 4 and the large P(j) for 
n — 5 is given by the monopole pairing interaction (to be discussed in Sec. 5) . One 
should be aware, however, that the large P(j)'s for an odd number of fermions are 
not necessarily related to P(0)'s of its even n neighbors, except for the connection 
given by the monopole pairing interaction. In many cases P(0)'s of systems with 
even numbers of fermions are much larger than P(j)'s of those with odd particle 
numbers. 

Figure 3 plots the P(0)'s for n = 4 up to j =33/2, for n = 6 up to j = 27/2, and 
the P(j)'s for n = 5 and n = 7 up to j = 27/2. It is seen that g.s. probabilities for 
even numbers of fermions in a single-j shell, and j g.s. probabilities for odd numbers 
of fermions in a single-j shell, stagger synchronously at an interval of 5j — 3 when 
j is small. When j is large, the P(0)'s and P(j)'s seem to saturate. 
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3.2 Bosons with spin / 

Figure 4 shows a few important P(/)'s versus I for four bosons with spin I [50]. One 
sees that the pattern of / g.s. probabihties for four bosons with spin I is very similar 
to that of four fermions in a single-j shell. For instance, the P(0)'s vs. I stagger at 
an interval of 6i — 3; the P(0)'s are dominant over other -P(/)'s except two I shells, 
1—2 or 8; and -P(/max)'s decrease with I. A new feature of / g.s. probabilities for 
four bosons with spin / is that the -P(0's are considerably large (small) when I is 
even (odd), i.e., the P(/)'s exhibit an odd-even staggering behavior. 

Next let us come to boson systems with spin / and particle number n. Here one 
should ask about the feature of P(/)'s for an odd number of particles: which angular 
momentum / g.s. dominance will appear, I — or I — 17 In systems with an odd 
number of fermions in a single-j shell, there are no 7=0 states, and one expects 
I — j g.s. probabihties to be large, as shown in Fig. 2. For bosons with spin I and 
odd n, however, one may have 7=0 states and thus it is not known a priori whether 
the I — or I — I g.s. dominance occurs in these systems. 

Figure 5 shows the P(/)'s and P(0)'s in boson systems with / —A and 6, and n 
running as large as possible. For odd n and / = 1, or for n = 3 and any odd / [48], 
there are no / = states; for odd n and / = 3, 5, 7 and 9, the / = states do 
not exist unless n > 15, 9, 7, and 5, respectively. For fifteen bosons with / = 3, 
nine bosons with / = 5, seven bosons with / = 7 and five bosons with / = 9 or 11, 
P(0) ~ 0% according to calculations by using 1000 sets of the TBRE Hamiltonian 
[50]. From Fig. 2 and these odd-Z cases it is concluded that the P(0)'s are usually 
much less than the corresponding P(0's when n is odd. On the other hand, the 
P(0)'s are mostly larger than the P(0's when n is an even number. 

These results indicate that the g.s. dominance is robust for systems with 
an even number of n, but not true generally if n is odd. In the latter case, it is 
observed that the g.s. dominance is easily lost. One therefore expects that the 
g.s. dominance is partly connected to an even number of particles. 

When one applies the TBRE Hamiltonian to sp bosons [51, 52], sd bosons [16, 
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17, 51, 53, 54, 55, 56, 57], and sdg [58] bosons (we shall come to these cases in 
subsections 4.6 and 5.2), the g.s. dominance is found also for an odd number of 
bosons. However, to a large extent, the dominant -P(O) therein is associated with 
the s boson condensation, which contributes around 40% to the ground state 
probability [53]. In other words, one should be aware that sp and sd systems are 
very special systems in which s boson condensation produces the g.s. dominance 
when the boson number is odd. Without s bosons, the P(0)'s of those systems with 
odd n would be drastically smaller and other I g.s. probabilities, such as those of 
I — I or I — /max, would be much larger. It would be very interesting to carry 
out systematic calculations of -P(/)'s for both even and odd numbers of bosons with 
various mixtures of spins. 

3.3 Many-j shells 

Let us exemplify the cases of fermions in many-j shells by using the pioneering 
works of Refs. [12, 13]. These authors calculated a few even-even nuclei in the sd 
shell (63 independent two-body matrix elements) and the pf shell (195 independent 
two-body matrix elements) with both neutrons and protons, with 1000 runs of a 
Hamiltonian by using the following ensembles: the TBRE, the random quasiparticle 
ensemble (RQE) which has an additional J dependence of the width [13], the RQE 
without monopole pairing (RQE-NP), and the RQE with splittings of single-particle 
energies (RQE-SPE). The results of four, six and eight neutrons in the sd shell, four 
protons and four neutrons in the sd-shell show that there are around 40% ~ 70% 
g.s. in these cases, although in these systems the / = states occupy a very small 
portion in the full shell model space. This result was a surprise to nuclear structure 
theorists, and has been attracting much attention since then. Table I shows g.s. 
probabilities by using different ensembles, and the percentage of / = states in the 
full shell model space. 

In Ref. [59], Horoi, Volya, and Zelevinsky extended the investigations of g.s. 
dominance to both even-even and odd-odd nuclei. They took random two-body 
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interactions uniformly distributed between —1 and 1. They checked the even-even 
nucleus ^'^Mg and the odd-odd nucleus ^^Al. They found that for ^^Mg the ground 
states are dominated by J = and T = (T is the total isospin of the state), 
and that states which satisfy the relation (— )^^"^ = 1 arc favored. For the odd-odd 
nucleus ^^Al the ground states are dominated by / = 1 and T = 0, and the favored 
states have {—Y'^'^ — — 1- In both cases the lowest T—0 is dominant. They also 
studied a simpler model, i.e., four nucleons in a single-j {j — 15/2) shell with both 
protons and neutrons. The situation is quite similar to that of ^^Mg: there is a 
similar predominance of the lowest /, T states, and a preference of the ground states 
with (— )^+^=l is easily noticed, while in the case of six protons and neutrons in 
a j = 9/2 shell the ground states with (— )^+^=— 1 are favored. Based on these 
examples, the authors concluded that the ground states are dominated by the states 
with / and T as low as possible, and quantum numbers satisfying a "selection" rule 
j-_-ji+T_^_^n/2 favored, where n is the number of particles. It is interesting to 
investigate more cases to see whether or not this observation is applicable for other 
systems. 

Now let us come to -P(/)'s of boson systems with spins more than one. The sd 
and sp boson systems were studied by Bijker and Prank [16, 17, 53, 54, 51, 55], and 
Kusnezov [52]. In Refs. [16, 17], (sd)" boson systems with n ranging from 3 to 16 
were considered by using random one- and two-body interactions, and a dominance 
of g.s. (~ 60%) was obtained for these cases. For sd bosons, the 7 = 2 (spin / 
of d boson) g.s probability and the /max g-s. probability arc also large, while other 
/ g.s. probabilities arc nearly zero; for sp bosons, P(/)'s arc very similar to the sd 
boson case except that 1 = 1 ground states (spin / of p boson) are favored instead 
of / = 2 ground states. 

3.4 Parity distribution in the ground states 

Another relevant quantity is parity distribution in the ground states in the presence 
of the TBRE Hamiltonian. For atomic nuclei, all the ground states of even-even 
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nuclei are observed to have positive parity, while those of odd-odd nuclei have both 
positive and negative parity with slightly more positive ones. The nuclei with odd 
mass numbers have almost equal number of cases with positive or negative parity. 
Table II presents the statistics of parity in the ground states of nuclei with mass 
numbers larger than 120. This table is based on the experimental data compiled in 
Ref. [60]. 

It is interesting to see whether a similar pattern appears in the presence of the 
TBRE Hamiltonian. A series of calculations was carried out in Ref. [61] for four 
model spaces: 

A Both protons and neutrons are in the f&pigg shell which corresponds to nuclei 
with both proton number Z and neutron number N ~ 40; 

B Protons in the fipigi shell and neutrons in the gid§. shell which corresponds 
to nuclei with Z ~ 40 and N ~ 50; 

C Both protons and neutrons are in the hiisids shell which corresponds to nuclei 

2 2 2 

with both Z and N ~ 82; 
D Protons in the qzds shell and neutrons in the hiisid^ shell which corresponds 

2 2 2 2 2 

to nuclei with Z ~ 50 and N ~ 82. 

These four model spaces do not correspond to a complete major shell but have been 
truncated in order to make the calculations feasible. These truncations are based 
on the subshell structure of the involved single-particle levels. It is noted that the 
number of states (denoted as D{I)) for positive and negative parity are very close 
to each other for all these examples. One thus expects that the probability of the 
ground states with positive parity is around 50%, if one assumes that each state of 
the full shell model space is equally probable in the ground state. 

We denote valence proton number and valence neutrons number by using Np and 
Nn, respectively. Because N and Z for closed shells are always even, nuclei with even 
values for both Np and Nn correspond to even-even type; nuclei with odd values for 
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{Np + Nn) correspond to odd-A type; and nuclei with odd values for both Np and 
Nn correspond to odd-odd type. As discussed above, the statistics for these three 
types are quite different. 

The calculated statistics of parity in the ground states obtained by using a TBRE 
Hamiltonian is given in Table III. For even-even nuclei (even values for both Np and 
Nn), it was noticed that positive parity is dominant in the ground states. For odd- 
A nuclei and doubly-odd nuclei, it is found that probabilities to have positive or 
negative parity in the ground states are almost equal with some exceptions. In 
general, there is no favoring for either positive parity or negative parity in the 
ground states of odd mass nuclei and doubly odd nuclei in the presence of random 
interactions. It is noted that these calculations are done for the beginning of the 
shell. For the end of the shell the results show a similar pattern. 

It was also found that the above regularities for parity distributions hold for 
very simple cases: single-closed two-j shells, one with positive parity and one with 
negative parity. The following shells have been checked: (2ji,2j2) = (9)7), (11,9), 
(13,9), (11,3), (13,5), (19,15), (7,5), (15,1). The statistics is very similar to the 
above results: The probability of ground states with positive parity is about 85% 
for an even number of nucleons, and about 50% for an odd number of nucleons. 

It is very interesting to note that for all even-even nuclei the -P(0+) is usually two 
orders of magnitude larger than P{0~). It would be very interesting to investigate 
the origin of the large difference in P(0) for positive and negative parity states, i.e., 
why the 0^ is not favored in the ground states. As is the case for an odd number 
of bosons with spin / (refer to Sec. 3.2), spin / = is not a sufficient condition to 
be favored in the ground states of a many-body system in the presence of random 
interactions, i.e., in order that the state is favored as a realistic ground state, not 
only 7 = but also positive parity is required. 

A simple system to study the parity distribution of the ground states in the 
presence of random interactions is the sp-boson system. First, it is noted that a 
sp-boson system with an odd number of particles n has the same number of states 
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with positive and negative parity; while for an even number of particles n there are 
slightly more states with positive parity (the difference is only n + 1). The parity 
for sp bosons is given by (— )^. The calculated results of Ref. [51] showed that when 
the number n of sp bosons is even, P(0) + P{n) ~ 99%, which leads to positive 
parity ground states dominance. When the number of n is odd, only about 50% of 
the ground states in the ensemble have / = 0, and about 50% have I — 1 or I — n. 
This leads to about equal percentages for positive and negative parity ground states. 
This pattern is very similar to that observed for fermion systems. 

3.5 Many-body systems interacting by a displaced TBRE 

While a TBRE is distributed symmetrically with respect to zero, an interesting 
question is what happens if one uses random interactions with only positive or neg- 
ative signs, or random interactions which are not distributed symmetrically around 
zero. This issue is both interesting and important because interactions in realistic 
systems, such as nuclei, atoms etc, are not symmetric around zero. Below we present 
results calculated by using a displaced TBRE. 

Let us firstly consider two arbitrary ensembles {Gj2-(jij2j3j4)} and {Gjr(jij2j3j4)}, 
which are related by a shift c: 

G'jTihhhk) = GjT{3i32h3i) + c, (17) 

where c is a constant. If {GjT{3i32j?,j4)} is symmetric aroimd zero and is described 
by Eq. (15), we call {GjT{3i323z3i)} the TBRE, and {Gjt (J1J2J3J4)} will be called 
a displaced TBRE with a displacement c. 

For fermions in a single-j shell, the results by using the ensemble {G'j} are 
exactly the same as those obtained by using {Gj}, except for a shift "^"~^^ c of the 
eigenvalue of the ground state. Therefore, a displacement of the TBRE is trivial in 
a single-j shell. 

In Ref. [62], the authors showed that for the /7/2P3/2 shell the displaced TBRE 
with an attractive average leads to the g.s. dominance and rotational spectra with 
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strongly enhanced B(E2) transitions for a certain class of model space. We shall 
discuss the results of Rcf. [62] in Sec. 6.2.3. 

In general, the role played by c in the displaced TBRE of Eq. (17) is very 
complicated [63]. For instance, both negative and positive displacements for an even 
number of fermions in many-j shells may favor (or quench) the g.s. probability, 
or produce a minor change. 

Recently, Johnson [64] also revisited the consequences of a displaced TBRE for 
^^Te and ^^Ca. He pointed out that one should be aware of a fact that the wave 
function obtained by the TBRE with an attractive average (i.e., c < 0) is quite 
close to that obtained by a negative constant value c for all the two-body matrix 
elements, if the width of TBRE Hamiltonian is 1 and c = —3 or —2. 

3.6 Effect of higher rank interactions 

As for the effect of interactions with rank higher than two, the only work was done 
by Bijker and Frank in Ref. [17]. These authors studied the case of sd bosons 
by adding three-body interactions to the one-body and two-body Hamiltonian (two 
one-body, seven two-body, and seventeen three-body interactions), where they used 
a scaling which depends on boson number n: For one-body terms this scaling is 
1/n; for two-body terms this scaling is l/(n(n — 1)); and for three body terms 
it is l/(n(n — l)(n — 2)). Their calculations showed that the inclusion of three- 
body interactions does not change the results in a significant way. When the boson 
number n is sufficiently large in comparison to the rank of interactions, the results 
are essentially similar. A study by using random ensembles of one- and two-body 
interactions also showed similar results to the case of pure three-body interactions. 
In other words, the basic features of sd boson systems do not change significantly 
due to the inclusion of three-body interactions. 



27 



3.7 Odd-even staggering of binding energies 

The odd-even staggering of binding energies is well known in nuclear physics. It is 
an evidence of the pairing interaction between like particles. Similar features were 
found in Ref. [13] for the angular momentum zero ground states by using the RQE 
Hamiltonian for four to ten neutrons in the pf shell. Johnson and collaborators 
made a least square fit of the binding energies, which was used to simulate Talmi's 
formula of binding energies within the framework of the generalized seniority scheme 
for atomic nuclei [14, 15], for the case of even numbers of neutrons in the presence 
of the RQE Hamiltonian. Then they calculated binding energies of systems with 
five and seven neutrons in the same shell by using the same sets of the RQE Hamil- 
tonian. Statistics of deviation of binding energies from Talmi's formula showed that 
deviations for an even number of neutrons are usually small and symmetric around 
zero, while those for an odd number of neutrons are large and positive. This is very 
similar to the situation of binding energies for atomic nuclei. 

In Ref. [18], Papenbrock, Kaplan and Bertsch considered the random Hamil- 
tonian for quantum dots or small metallic grains which conserve total spin. The 
results of Ref. [18] showed that even a purely random two-body Hamiltonian can 
give rise to the odd-even staggering of binding energies. 

3.8 Summary of this Section 

In this Section we showed the robustness of the g.s. dominance for the case of 
even numbers of fermions in a single- j or many-j shells, and also for the cases of 
(both even and odd numbers of) sd or sp bosons. There are only very few counter 
examples for the case of four fermions in a single- j shell (j = 7/2 and 13/2). For 
neutron-proton systems, it was found that the P(0) is not dominant for odd-odd 
systems [59, 61]. For bosons with spin /, there are many cases in which the P(0)'s 
are very small in comparison to other P(J)'s when the boson number n is odd, 
which suggests that the g.s. dominance is partly connected to an even number of 
particles. In other words, the g.s. dominance in many-body systems is a robust 
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feature associated with even numbers of particles. 

The parity distribution in the ground states calculated by using the TBRE 
Hamiltonian is found to be robust: for the case of even numbers of both valence 
protons and neutrons the positive parity is always dominant (around 80%), while 
positive parity and negative parity are almost equally probable for other cases. For 
atomic nuclei, all even-even nuclei have positive parity g.s. and the other cases have 
either positive or negative parity g.s. with almost equal probabilities (refer to Table 
II). 

A displaced TBRE Hamiltonian may give a very different pattern of -P(/)'s 
for the case of fermions in many-j shells in comparison to those obtained by a 
TBRE Hamiltonian. This means that the / g.s. distribution for a displaced TBRE 
Hamiltonian is a much more complicate issue. 

The effect of the rank of the Hamiltonian was studied in Ref. [17] for srf-boson 
systems. The calculated results of Ref. [17] showed that basic features obtained 
by using the random Hamiltonian do not change significantly even if one includes 
three-body interactions. 

Systematic odd-even binding energy differences were first discussed in Ref. [13] 
for a few nucleons in a sd shell, and revisited in Ref. [18] for finite metallic clusters 
by using the TBRE Hamiltonian. These calculations concluded that an odd-even 
staggering arises from the TBRE Hamiltonian: stronger binding energies for systems 
with even numbers of particles are typically obtained in numerical simulations. How- 
ever, the g.s. dominance and odd-even staggering are not necessarily two facets of 
the same thing. 
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4 Spin / ground state probabilities of simple sys- 
tems 



In this Section we discuss P(J)'s of simple systems. Here the appellation "simple 
systems" means either the eigenvalues of the systems are linear in terms of two-body 
matrix elements or one can classify the ranges of two-body matrix elements (because 
of some specific features of the systems) in a simple way (see Sec. 4.4). 

In this Section we shall first come to a few systems in which the eigenvalues are 
linear combinations of two-body matrix elements. One feature for a state to have 
a large g.s. probability, and an empirical approach to predict -P(/)'s which were 
introduced in Refs. [47, 56, 57, 58], will be discussed. For three and four fermions 
in a single-j shell (j < 7/2) and d boson systems, one can calculate -P(/)'s exactly 
by using the geometry of eigenvalues in the two-body matrix elements. This recipe 
was introduced in Ref. [65], and will be discussed in this Section. The evaluation of 
P(/)'s of sp and sd bosons by a mean-field approach introduced in Refs. [51, 53, 54] 
will be also discussed. 

Due to the simplicity and specific properties of the systems discussed in this 
Section, P(/)'s of these systems can be predicted very well. Thus one can say that 
P(/)'s of these systems in the presence of random interactions are satisfactorily 
described. 

4.1 Fermions in a single-j shell with j < 7/2 

The eigenvalues of states of fermions in a single-j shell with J < | can be written 
in terms of linear combinations of the two-body matrix elements. Here we discuss 
only the case of four fermions in a j = | shell. Discussions of systems with n = 3 
and j < I can be found in Ref. [58] . 

For a j — 7/2 shell with four fermions, all the states are labeled by their total 
angular momenta / and their seniority quantum numbers (v). The eigenvalues £'/(„) 
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are as follows (see Ref. [66]): 
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In Eq. (18), bold font is used for the largest and italic for the smallest amplitudes 
in an expansion in terms of Gj. Eq. (18) can be rewritten as follows, 

Em^T.^mGj, (19) 
J 

where j3 representing all the necessary additional quantum numbers to label the 
state. By using the TBRE Hamiltonian described by Eq. (15) and the eigenvalues 
given by Eq. (18), it is easy to obtain the probability, P{I), for each / ground state. 
/ g.s. probabilities for four fermions in a j = 7/2 shell are shown in the row "TBRE" 
of Table IV, and are obtained by 1000 runs of the TBRE Hamiltonian. 

One can also predict the / g.s. probabihty without running the TBRE Hamilto- 
nian. For example, the P{0) is determined by 

fdGofdG2fdGjdG,p{Go)p{G2)p{G^)p{G,)\. x, (20) 

where I{v) 7^ 0(0), and the subscript "(Z^j Q!o(o)Gj < J2j c>ii{v)Gjy^ is the require- 
ment for Gj's which take values from —00 to 00. 

The row "predl" of Table IV corresponds to probabilities predicted by an integral 
for each 7+ state similar to Eq. (20) for the 0+ state. Probabilities calculated by 
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using the TBRE and those predicted by using integrals hke Eq. (20) are consistent 
within statistical fluctuations. 

One easily sees from the present example that a state with one or more largest (or 
smallest) a/^^,-), for which we used bold (or italic) font in Eq. (18), has a very large 
probability to be the ground state (or the highest state). The quantum numbers 
I{v) of these states are: i"(v)=0(0), 2(4), 4(4), 8(4). P(/)'s of states without the 
largest and/or the smallest aj for a given J are very small. As we shall see from 
the discussions of next subsection, this regularity is the feature for a certain state 
to have a large probability to be the ground state. 

The regularity that large -P(/)'s are related to the largest or the smallest co- 
efficients among the fixed two-body matrix elements, discussed above, is actually 
very general. Now let us look at linear combinations of the random numbers which 
follow the Gaussian distribution with the average being zero and the width being 
one. Suppose that F{k) be a set of linear combinations of Gj: 

F{k)=Y.aiGj, k = l,2,---IC, (21) 
J 

where k is used to specify each F, and K. is the total number of different F's. One 
can prove that distribution functions of random F{k) are 

''''''' -vkr^i-'-W)' 

If in Eq. (21) is the largest (or the smallest) among all the a{ {k = 1, ■ ■ -/C), 
the probability of F{m) being either the smallest or the largest number is large. To 
show this, let us look at 

J^{k) = F{k) - F{m) = (ai' - a^) Gj, + i^Y.^ {ai - ai) , (23) 

where k ^ m. The right hand side of Eq. (23) has two terms, both of which are ran- 
dom numbers which follow the Gaussian distribution. The value of ^a^' — a^^ Gj/ 
is negative or positive for all ^(A;)'s, and thus effectively produces a "simultaneous" 
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shift to either negative side or positive side for all ^(/c)'s, depending on the sign 
of {a'j! — Gji. Therefore, all of the functions J-'(kys have large probabilities to 
be both negative and positive, i.e., F{m) has a large probability to be either the 
smallest or the largest. 

If there are two or more coefficients a'^ (J' = 0, 2, ■ ■ ■ 2j — 1) which are the largest 
or smallest for different functions F{k), the probability of finding F{m) as either 
the smallest or the largest is expected to increase. 

By using Eq. (18) and Eq. (22), the distribution width, gi ^, for each state of 
four fermions in a j = 7/2 shell are fisted in the last row of Table IV. It is easily 
noticed that there is no correlation between the P{I) values and their corresponding 
g'j^^f. For example, P(0) (~ 20%) is much larger than P(5) (= 0) and P(6) (~ 0) 
although 5'7(^)=o(o) is much smaller than gi(v)=b{i) and gi{v)=&{2) ■ 

4.2 An empirical approach 

Because P{I) 's in Eq. (20) cannot yet be determined by a simple analytic procedure, 
one has to evaluate this integral numerically. It is therefore desirable to find a simple 
alternative method to evaluate / g.s. probabilities. 

In Refs. [56, 58], such a substitute was given. The idea is based on the ob- 
servation discussed above: The state with a sizable / g.s. probability involves the 
largest and/or smallest Q;/(^-) with J fixed. Thus the / g.s. probability might be 
proportional to the number of the largest and/or smallest aj(^.^y Let M[ be the sum 
of numbers of the smallest and the largest q;/^^^ with a fixed J for a certain /. Then 
the / g.s. probability is approximately given by 

P{I)=N'i/Nm, (24) 

where Nm — '^N — 1 with N the number of two-body matrix elements. Note that 
= 2N — 1 is used instead of 2A'", because all aj^'^{I ^ 0)'s are (there is no 

smaUest a/(^°), and that Y^jN'j = N^- 
^We shall define another width, aj, in Sec. 5.6. 
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Now we exemplify this empirical approach by using four fermions in a j = 7/2 
shell Here iV = j + i = 4, and iV„ = 2 x 4 -1 = 7. From Eq. (18), it is easy to find 
the largest (or the smallest) a/^^-j for different I{v) states but fixed Gj: 

3 

for Go Q^o(o) = — is the largest ; 
2 



7 10 

for G2 q;4(4) (or al(^^^) = — (— ) is the largest (smallest) ; 

4 42 

for G4 a2(4) (o^ ci;4(4)) — ( ^ ) largest (smallest) ; 

12^ 
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for Ge q;8(4) (or 0:2(4)) — (77) is the largest (smallest) . 



We thus find that 



Aq = 1 given by the largest q;2(o); 

J\f2 — 2 given by the largest 0:2(4) smallest 0:2(4) j 

A/'4 = 2 given by the largest 04(4) and the smallest a^^^) ; 

= 2 given by the largest 08(4) and the smallest 0:5(4) . 

According to the above empirical approach, we therefore predict that P(0) is 1/7, 
and P{I) with 7 = 2, 4, 8 and v = 4 is 2/7. These predicted / g.s. probabilities are 
given in the row "pred2" of Table IV. A very reasonable agreement is easily noticed 
between the results obtained by running the TBRE Hamiltonian, those obtained by 
multiple integrals such as Eq. (20), the solutions by using geometry method which 
will be discussed later, and the predicted values by the present empirical approach. 

The largest (smallest) o;/(„)'s with fixed J correspond to the ground (highest) 
state when Gj — —1 and other Gji^jS are zero. Below the ground (highest) state 
for one of Gj = — 1 and other Gji^j — will be discussed instead of the largest 
(smallest) oi^^^''s, for the sake of convenience. 

Now let us study another case: d boson systems. Similar to fermions in a single-j 
shell for small j (j = | or j = |), the relation between the two-body matrix elements 
and the eigenvalues for d-boson systems is also linear. The two-body Hamiltonian 



34 



of a d-boson system is given by 

Ha=j:lc,{sd^y^\ddtK (25) 

From Eq. (2.79) of Ref. [32], one obtains 

E^E^ + a'^ridina - I) + P' K(nrf + 3) - v{v + 3)] + 7' [/(/ + 1) - 6nJ , (26) 

where £^0 contributes only to binding energies, and Ud is the number of d bosons. 
Eq. (26) can be rewritten as follows, 

E{v, TiA, I) = E'oiTid) - P'v{v + 3) + 77(/ + 1). (27) 

The a', P' and 7' in Eq. (26) are linear combinations of Cq, C2 and C4. Prom Eq. (2.82) 
of Ref. [32] , one obtains 

a' = ^(4c2 + 3C4), 

/3' = ^(7co-10c2 + 3c4), 

7' = ^(-C2 + C4). (28) 

Substituting these coefficients 7' into Eq. (27), and taking the two-body 
matrix elements co,C2 and C4 to be the TBRE defined by Eq. (15), one easily 
calculates / g.s. probabilities which arc shown in Fig. 6. It is easy to notice that 

1. The P(/max)'s are almost constant (around 40%) for all (< 4); 

2. The P(0)'s and P(2)'s arc periodic, with a period 6{nd)=Q. 

3. All the P(/mc«), P(0) and P(2) are near to 0, 20%, 40%, or 60%. Other 
P(/)'s are always zero. 

Now let us explain these observations by using the empirical approach introduced 
in Refs. [56, 58]. Prom Eq. (27) and Eq. (28), 

Co = -1, C2 = C4 = : E{v, Ua, I) = -£^o(^d) + + 2)' 

1 1 

C2 = -1, Co = C4 = : E{v, UA, I) = £^oM - j'^i'v + 3) + j^Hl + 1); 

C4 = -1, Co = C2 = : E{v, UA, I) = E'oiud) + ^v{v + 3) - ^/(/ + 1), (29) 
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where EQ^na) is a constant for all states. Based on Eq. (29), one obtains Table V, 
which presents the angular momenta giving the largest (smallest) eigenvalues when 
Cl = —1 {L =0, 2, 4) and other parameters are for d boson systems. In Table V, 
AC is a natural number, and > 3. These angular momenta appear periodically, 
originating from the reduction rule of U(5)— s>S0(3). 

One can use Table V to predict the P{I) for d-boson systems. For example, 
jVp = 3 and A/}^^ = 2 forn = 6A;. According to the empirical approach of Eq. (24) 
^ , the predicted P{I) = Afj/N^, where = 5. Thus the predicted P(0) = 60% 
and -P(/max) = 40%, and all other -P(/)'s are predicted to be zero for n — 6k. 

A comparison between the P{I) values predicted by Table V and those in Fig. 
6 shows that a certain P{I) is large if one state with angular momentum I involves 
the largest and/or the smallest (Eq. (19)) for a given L {L = 0,2,4). 

For (i-boson systems one also finds that g.s. probabilities are very close to zero 
periodically when n = 6k ± 1 (k is a natural number). These counter examples of 
the g.s. dominance are also predicted by Table V: if n = 6k± 1, the / = states do 
not produce the largest and/or smallest eigenvalues when one of the parameters q 
is —1 and others cl' {V ^ L) are zero. This feature is consistent with the discussion 
of Sec. 3.2, where it was seen that the P(0) for an odd number of bosons with spin 
I is usually not dominant. 

4.3 Predicted P{iys based on geometry of the eigenvalues 

Chau, Frank, Smirnova, and Isacker [65] showed that the simple systems discussed 
in last subsection can be projected to a polyhedron on the axes defined by two-body 
matrix elements. These authors related the above ct^^ with the largest possible 
convex in W dimensions, where 1^ = — 1. In this subsection we shall discuss 
predicted -P(/)'s based on geometry of the eigenvalues, which was suggested in Ref. 

m- 

*Note that when one searches for the smallest eigenvalue with co = — 1 and C2 = C4 = in 

case A of Eq. (29), one finds that many states with different / arc degenerate at the lowest value. 
Therefore, again, = 2A^ — 1 = 5 is used in predicting P(/)'s by the formula P{I) = Afj/Nm- 
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Chau et al. discussed d bosons and four fermions in a j = 7/2 shell. Here let us 

take five d bosons as an example, where 1^ = 3 — 1 = 2. The procedure suggested 
by Chau et al. is as follows, 

1) One rewrites the eigenvalues of Eq. (26) as 

£ns,v,I = Ci + '^Cni,v,l{cL — (30) 
L 

where £ns,vj = ^'na,^,/^;^ = ^K,,^,/, is called the "scaled energy" [65]. 

2) The scaled energy of an arbitrary eigenstate is represented as a point in a 
plane spanned by (c2 — Cq) and (C4 — Cq). All points (corresponding eigenvalues) 
were found to be confined to a compact region with the size of one unit in each 
direction. In the case of d bosons, the / g.s. probability is related to each angle at 
the corresponding vertex i by 

p(/). = i-|. (31) 

where p{I)i is the probability of / to be the ground state contributed from the 
vertex i. The eigenstates for which corresponding points are not vertices can not 
be the ground states for the TBRE Hamiltonian. In Fig. 7 the twelve solid circles 
have one-to-one correspondence to the twelve states with different / and v for five 
d bosons. The angles (i = 1, • • • 4) at vertex are also labeled in Fig. 7. 

3) One finally sums p{I)i over all the i vertices to get the total P{I). 

This method is appropriate to discuss -P(/)'s of systems in which the number of 
two-body matrix elements is not large (3 or 4) (where the angle is relatively easy to 
evaluate), and meanwhile the Hamiltonian is diagonal. For more complicated cases 
this approach should be generahzed. Similar to Eq. (20), the approach of Ref. [65] 
presents exact -P(/)'s without using random interactions. 

4.4 Mean-field method 

P(/)'s of sp bosons in the presence of the TBRE Hamiltonian were first addressed 
by Kusnezov in Ref. [52]. He addressed -P(/)'s by using random polynomials. The 



37 



procedure is as follows: First, choose sp a system with a large boson number for 
which the dimension of the Hilbert space is large. The Hamiltonian is tlien reduced 
to a tri-diagonal form with the trial Lanczos state for which the number of p bosons 
is equal to I. Next one constructs the off-diagonal and diagonal matrix elements in 
terms of Up/n and I / Ima^- 

In the limit of large boson number n, the lowest eigenvalues for each / states 
can be found in terms of these off-diagonal and diagonal matrix elements of the 
tri-diagonal matrix, and they can furthermore be written in the form of a parabolic 
function of Up/n, with the coefficients determined by matrix elements of interactions 
and ///max- Then one is able to analyze the properties of each term to evaluate / 
g.s. probabilities. The advantage of the approach in Ref. [52] is that it avoids 
the diagonalization of matrices. In general, the lowest eigenvalues are not quadratic 
polynomials, but have a more complicated form and should be studied more carefully, 
as pointed out in Ref. [54]. The results of Ref. [52] is consistent with the mean-field 
approach of Refs. [53, 51, 55], which will be introduced in this subsection. 

In Refs. [51, 53], Bijker and Frank suggested a mean-field analysis to predict / 
g.s. probabihties in the vibron model and the sd IBM. They used the connection 
between potential energy surfaces of the Hamiltonian and geometric shapes. Let us 
discuss below only the case of sp bosons with a schematic Hamiltonian, because the 
philosophy for a general vibron Hamiltonian and the sd IBM Hamiltonian is similar. 

The schematic vibron Hamiltonian of Ref. [51] is given by Eq. (13) in Sec. 2.1. 
The range of x in Eq. (13) is from — | to The coherent state of the vibron model 
is given by 



The equilibrium shape is obtained by calculating one- and two-order derivatives 




(32) 



where a e 0) f • The potential surface is given by the expectation value of the 
vibron Hamiltonian in the coherent state, i.e., 

E{a) — -sinx sin^cc -|- cosx sin^o;. (33) 
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of E{a) with respect to a, and the results can be classified into the following three 
classes: 



(i) ao = 0, -- < X < 2 ; 
(n) cos2Qio = cotx, ^ < X < 

(m) ao = -,— <x<Y- ^ ^ 

The case (i), which corresponds to spherical symmetry and produces only 
ground states, occupies a portion of (37r/4)/(27r) and contributes 3/8=37.5% to 
P(0); the case (n), which gives ground states when f < x ^ f and /max = n 
ground states when f < x < T' contributes (l7r/4)/(27r) = 12.5% to both P(0) 
and P{n); the last case (iii), which produces g.s. (if n is an even number) or 1 
g.s. (if n is an odd number) when vr < x ^ and produces /max = n g.s. when 
^ < X ^ TT, therefore the case (iii) contributes ^/(27r)=25% to -P(O) for even n 
or 25% to P(l) for odd n, and contributes |/(27r)=12.5% to P{n). To sum over 
these three cases, one obtains P(0) = 75% and P(n) = 25% for even n and that 
P(0) = 50%, P(l) = 25%, and P(n) = 25% for odd n. 

P(/)'s obtained by using the TBRE Hamiltonian for vibrons, shown in Fig. 8, 
are very close to this simple prediction. A detailed mean-field analysis for sp bosons 
can be found in Refs. [51, 53, 54]. 

This mean-field approach was also applied to predict P(/)'s of sd boson systems 
[53, 55]. Recently, Kota applied this approach [67] to analyze probabilities of dif- 
ferent irreducible representations in the ground states obtained by using the TBRE 
Hamiltonian of the IBM. 

4.5 Summary of this Section 

In this Section we first showed that for systems in which eigenvalues are linear 
combinations of two-body matrix elements, one can apply an empirical approach 
suggested in Refs. [47, 56, 58] or an approach based on geometry of eigenvalues 
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suggested in Ref. [65] in order to predict -P(/)'s without using random interactions. 

The examples include d boson systems, three and four fermions in a j=5/2 or 7/2 
shell. 

We also showed that one can explain P(/)'s of systems (such as sp- and sci-boson 
systems) in which one can classify two-body matrix elements in a simple way. P{I)^s 
for sp and sd bosons in the presence of the TBRE Hamiltonian were evaluated in 
Refs. [54, 51] based on a mean-field method. The g.s. dominance for sp bosons 
were also discussed in Ref. [52] by Kusnezov based on random polynomials. The 
results of these two approaches are consistent, but the mean-field approach has a 
more transparent picture. 

On the other hand, one is unable to predict P(/)'s of the systems discussed in 
this Section without a priori knowing the relations between the eigenvalues and 
the two-body matrix elements or the specific features of their eigenvalues. In this 
sense, the behavior of P(/)'s for such examples discussed in this Section has not 
been understood at a fundamental level. The mean-field approach for the sp and sd 
bosons [54, 51] has not been successfully generalized to other cases. 

Because the problem of the g.s. dominance is very difficult, one may leave 
this embarrassing situation to the future and proceed by describing i-*(/)'s for more 
complicated systems in which the eigenvalues are not linear combinations of two- 
body matrix elements and one is unable to classify the two-body matrix elements in 
the way discussed in Sec. 4.4. Fortunately, the empirical approach of Sec. 4.2 was 
found to be applicable to these cases, after very slight modifications. This will be 
discussed in the next Section. 
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5 Spin / ground state probabilities of complicated 

systems 



In this Section we shall discuss P(/)'s of compUcated systems. Here the appellation 
"compUcated systems" means the eigenvahics of the systems are not hnear in terms 
of two-body matrix elements and one cannot classify two-body matrix elements in 
a simple way (as done for sp or sd bosons in last Section). 

Although there have been no simple understandings of the g.s. dominance in 
complicated systems so far, the empirical approach discussed in last Section was 
found to be reasonably applicable to prediction of -P(/)'s, after slight modifications. 
In this Section we shall first discuss results along this line. 

We shall first consider fermions in a single-j {j > 7/2) shell, next go to fermions 
in many-j shells and boson systems. An argument of the empirical approach will 
be discussed schematically. Then we shall discuss P(/max) 's which were found to be 
considerably large for fermions in a single-j shell and boson systems. 

In this Section we shall also discuss an argument of the g.s. dominance for 
four fermions in a single-j shell or four bosons with spin / [68], studies of the 
g.s. dominance based on the distribution width of the eigenvalues for states with 
different /, and those based on time reversal invariance of the Hamiltonian, large 
overlaps between g.s. wavefunctions of systems with mass number differing by 
two, results of ground states for spin-1/2 fermions, etc. 

5.1 Fermions in a single- j shell 

Let us take the case of four fermions in a single-j shell again. This case was first 
studied by Mulhall, Volya and Zelevinsky [45, 46] by using two-body interactions 
which distribute uniformly between —1 and 1. However, ground state probabilities 
discussed by these authors are not -P(/)'s, but actually those of centroids of states 
with given I (denoted as V{I) in Sec. 6). For instance, they predicted 50% of 
-P(-^min) and another 50% of -P(/max), which is far from the observation but is the 
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behavior of energy centroids with spin /. 

In Rcfs. [56, 58] the empirical approach was gcnerahzed as follows: First, one sets 
one of the two-body matrix elements Gj to —1 and all others to zero. Then one finds 
which angular momentum / gives the lowest eigenvalue among all the eigenvalues 
of the shell model diagonalization. Suppose that the number of independent two- 
body matrix elements is N , then the above procedure is repeated N times. Next, 
among the N runs one counts how many times (denoted as Hi) a certain angular 
momentum I gives the lowest eigenvalue among all the possible eigenvalues. Finally, 
we predict that the probability of / g.s. is given by 

P{I)^Ni/N. (35) 

In Eq. (24) of Sec. 4.2 we used A/} which is the number of both the largest and the 
smallest eigenvalues for Gj — —Sjj/, and which is equal to (2A^ — 1) for four 
fermions in a j = 7/2 shell ^, in order to predict the values of P{I) empirically. 

Below let us use Eq. (35) to predict -P(/)'s of comphcated systems, namely, we 
do not use the largest eigenvalues for Gj — —5jji. The reason of this modification 
is that the largest eigenvalues are usually (exactly or nearly) zero for matrices corre- 
sponding many spin /'s, especially for many-j shells or a single-j shell with large j. 
To have the empirical approach as simple as possible, we shall use only the lowest 
eigenvalues with one of the Gj's being —1 and the others being switched off, for 
the case of fermion systems in a single-j shell for large j, many-j shells, sd- and 
sdg-hoson systems. 

For fermions in a single-j shell, tables of the angular momenta / which give the 

lowest eigenvalues for n — 3-7 were presented in Refs. [56, 58]. The cases of n = 4 

is given in Table VI as an example. Here we also mention two systematic features. 

The first concerns the quadrupole matrix element (6*2) term. It has been well known 

that the monopole pairing interaction always gives an 7 = ground state for an even 

^One can also use Eq. (35) for four fermions in a j = 7/2 shell. The predicted -P(O) = f (2) = 
P(4) = P(8) = 25% if one takes only the largest a's. These predicted -P(/)'s are also reasonably 
consistent with other results. 
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number of fermions in a single-j shell, and em I = j ground state for an odd number 
of fermions in a singlc-j shell when Go is set to be —1 and others 0. However, little 
was known about the role of the G2 matrix element in a single-j shell. It was found 
in Ref. [58] that the J = 2 pairing interaction always gives an / = n ground state 
for an even number of fermions, and an / = j — (n — l)/2 ground state for an odd 
number of fermions, when G2 is set to be —1 and others 0. Another regularity is 
that interactions Gj — —1 (all others are zero) with J — 2j — 3,2j — 5 ■ • ■ {j > 9/2) 
produce / = /max — 8, /max " 16 • • • g.s. for n — A, Gj — —1 with J — 2j — 3 
{j > 13/2) produces an / = /max — 12 g.s. for n — 5, and Gj — —1 with J — 2j — 3 
{j > 15/2) produces an / = /max — 20 g.s. for n — 6. No understanding of these 
features are available. 

Now we exemplify the applications of Table VI for four fermions in a single-j 
shell by the case of j = 9/2, where iV = j + l/2 = 5. Here one finds according to 
the second row of Table VI that A/q = 3, and A/4 = A/^^^ = 1. The predicted values 
of P(/)'s are P(0) = 60% and P(4) = P(/max) = 20% while all other predicted 
P(/)'s are zero. For 1000 runs of the TBRE Hamiltonian [56], the P(0) = 66.4%, 
P(4) = 11.8% and P(/max) = 17.9% while all other P(/)'s are close to zero. The 
agreement between the predicted P(/)'s by the empirical method of Refs. [56, 58] 
and those obtained by diagonahzing the TBRE Hamiltonian is thus very good. 

This empirical method can be also applied to predict P(/)'s of odd numbers of 
fermion systems. For example, for five fermions in a j = 9/2 shell, one finds from 
the first row of Table IV in Ref. [58] that Afj = 2, Afs/2 = A/^/a = A/}^,, = 1. 
The predicted values of P(/)'s are thus P(3/2) = P(5/2) = P(/max) = 20% and 
P(j) = 40% while aU other P(/)'s are predicted to be zero. P(/)'s obtained by 1000 
runs of diagonahzing the TBRE Hamiltonian are: P(3/2) = 20.5%, P(5/2) = 15.5%, 
-P(/max) = 18.4%, and P{j) = 33.9% while all other P(/)'s are close to zero. Good 
agreement is also obtained between the predicted values of P(/) 's and those obtained 
by using the TBRE Hamiltonian. 

Figure 9 gives a comparison between the predicted P(0)'s (open squares), which 
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are obtained by using Table VI of this paper and Table V of Ref . [58] , and those 
obtained by diagonalizing the TBRE Hamiltonian (solid squares) for n = 4 and 6. 
It can be seen that a good agreement is obtained for fermions in a single-j shell 
for both small and large j. The predicted g.s. probabilities exhibit a similar 
staggering as those obtained by diagonalizing the TBRE Hamiltonian. 

It is interesting to note that the P(0)'s of four and six fermions can also be fitted 
by empirical formulas: 



forn = 4:P(0) = K^^±ll/^xl00%, 



2 



1 if 2j = 3m 

if 2j + 1 = 3m ; 

-1 if 2j - 1 = 3m 



fo,„ = 6:P(0) = aM±lVl^xl00%, (36) 

^ 2 

where the " [ ] " means to take the largest integer not exceeding the value inside. 
These empirical formulas are interesting because they present scenarios for very 
large-j cases where it would be not possible to diagonahze the TBRE Hamiltonian. 
A comparison between the values predicted by these formulae and those obtained 
by diagonalizing the TBRE Hamiltonian is also given in Fig. 9. 

A correlation can be found in this figure between P{j) for an odd number of 
fermions in a single-j shell and -P(O) for a neighboring even number of fermions 
in a single- J shell. This correlation originates from an attractive monopole pairing 
interaction: Go = — 1 gives the g.s. for fermions with an even number n of particles 
and j g.s. for those with the neighboring odd n. Namely, at a large probability, 
random interactions with attractive Go give g.s. for the case of even n and j g.s. 
for the case of odd n. If one switches off Go term, there will be in general much 
fewer cases for which the same sets of random interactions give g.s. for the case 
of even n and j g.s. for the case of the neighboring odd n. 

5.2 Fermions in many-j shells, sd- and sc^^'-boson systems 

Now we discuss P(/)'s of fermions in many-j shells and bosons with many spin Ts 
by using the empirical approach of Eq. (35). This formula remains the same here 



44 



except that the Gj should be replaced by the general two-body matrix elements 

GjTijlj2,jd4)- 

Let us first exemplify a two-j (j = |, |) shell with n=4 to 7. A comparison of 
the predicted P(/)'s by using the empirical formula of Eq. (35) with those obtained 
by diagonalizing the TBRE Hamiltonian of fermions in a two-j (j = |, |) shell with 
n—4 to 7 is shown in Fig. 10. One sees a reasonable agreement ^. 

Now let us study sd-boson systems similarly. Table VII presents spins / which 
appear in the lowest states when one of the above sd-boson parameters is set to be —1 
and others 0. We predict, according to Table VII and the procedure of the empirical 
method, that only 7 = 0, 2, and 2n g.s. probabilities are sizable, while other I g.s. 
probabilities are close to zero. Fig. 11 shows a comparison of the predicted P(/)'s 
and those obtained by diagonalizing the TBRE Hamiltonian of srf-boson systems, 
with boson numbers ranging from 6 to 16. One can sec a reasonable agreement 
between P(/)'s obtained by using the TBRE Hamiltonian and those predicted by 
the empirical formula of Eq. (35). 

Because the empirical formula of Eq. (35) reasonably predicts P(/)'s of both 
simple and complicated systems, it is also called an empirical rule of P(/)'s in this 
paper. 

5.3 Simple argument of the g.s. dominance for four fermions 
in a single-j shell and four bosons with spin / 

A simple argument of g.s. dominance for four fermions in a single-j shell and four 

bosons with spin / was recently discussed in Ref. [50]. The essential point is that 

^For fermions in many-j shells, the number of two-body matrix elements is usually large. In 
such cases, especially for odd-fermion systems, there are "quasi-degeneracy" problems in count- 
ing A/}: sometimes the lowest cigcnvahic is quite close to the second lowest one when one uses 
Gj{jij2,j3j4) = — 1 and others 0. For such two-body matrix elements, one should actually intro- 
duce an additional "rule" in order to have a more reliable prediction. Namely, it is not appropriate 
to count J\fi in the most naive way. In order to avoid confusions, however, we do not modify 
the way to count A// in such cases throughout this article. It is noted that / = | in Fig. 10(b) 
and / = I in Fig. 5(d) are cases of "quasi-degeneracy". Improvement of agreement between the 
predicted -P(/)'s and those obtained by diagonalizing the TBRE Hamiltonian can be achieved by 
appropriately considering the above "quasi-degeneracy" . 
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there is only one / = state with non-zero eigenvalue among many 1 = states of 
these systems, while those of / 7^ are scattered with smaller values in magnitude 
when only one of the two-body matrix elements is taken to be —1 and others are 
zero. This feature can be proved by constructing the / = states by using a pair 
basis [68]. We shall discuss below for the case of four fermions in a single-j shell. 
According to the empirical rule, the / = states have a large probability to be the 
ground states for the TBRE Hamiltonian. 

The single-j Hamiltonian for which one of the Gj's is —1 and others are zero is 
defined by 

^ = ^. = -E^M^^M- (37) 

M 

We define 

: /)) = (A('^i)t>l('^2)t)^'^ |0) , (38) 



where |)) means that the basis is normalized, A*^*"*)^ is defined by Eq. (3) without 
isospin degree of freedom, and Nr-^r^i is given by 



Nr,r2l = 1 + <^rir-2 " 4(2ri + l)(2r2 + 1) < 



The matrix elements of Hj are as follows (see Ref. [69]), 




(39) 



{{ry^ : I\Hj\r,r2 : I)) = ^==^== ^ Ur.r^JRUr'.r'^jR , (40) 



with 

Urir2JR — ^Jri^Rr2 + {~y^Jr2^Rri — 4:LiL2JR < 




(41) 



where fi is a short hand notation of \^2ri + 1. 

Let = \JJ : 0)), \^k) = \KK : 0)) - {{KK : 0|JJ : 0))|JJ : 0)) {K ^ J). 
The new basis constructed by |$ft-)'s are orthogonahzed with respect to only 
and not normalized except Similar to Eqs. (41) and (39), one has 

{{KK : 0\Hj\K'K' : 0)) = {{KK ■.0\JJ : 0)) {{K'K' ■.Q\JJ : 0))Njjo, (42) 
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where K' = 0,2, ■ ■ ■ ,2j - 1. 

By using this formula, one easily confirms that all matrix elements {^x\Hj\^k') 
are zero except when K = K' = J. Namely, {{J J : 0\Hj\JJ : 0)) = Njjo = Ej^-'^ is 
the only non-zero eigenvalue for the / = states of four fermions interacting by Hj 
(This procedure is also applicable for four bosons with spin /). From the sum rule of 
diagonal matrix elements [71] one obtains that X]j -E'/*^'^'' — —^n{n—l)Dj^^ — —6Dj^\ 
where n is the number of fermions in the system (n — 4 here), and D^^hs the number 
of states with spin I. For / 7^ states the eigenvalues for Gj — —Sjji are scattered 
while for 7 = there is only one non-zero eigenvalue which is on average larger in 
magnitude than eigenvalues for states with other 7's except for very small j ^. This 
suggests that A/q is larger than A/} (/ 7^ 0). According to the empirical rule, the 
g.s. probability is larger than that of all other / states. 

One sees that the increase of P(0) "coincides" with that the number of J = 

states for four fermions in a single-j shell or that for (denoted as -Do"*) four bosons 

with spin I, respectively. Dq^ and Dq^ take their values 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 

4, • • • for 2j=3, 5, 7, 9, 11, • • • or 0, 1, 2, 3, 4, • • •, etc. Although it is still difficult 

to prove this correlation, below we give an argument based on the empirical method 

of Refs. [56, 58]. Taking a system of four fermions as an example, the relation 

J2j Ej^^ — —&D^p means that the magnitude of the sum of -Bq increases with 

Dq \ suggesting that Nq increases with D^"^ simultaneously. Because P(0) = Mq/N, 

a regular increase of D^^ of four fermions [49, 61] produces a regular staggering of 

the P(0)'s. It is noted that a similar "coincidence" of the staggering of the P(0)'s 

with increasing number of /=0 states was observed in Refs. [46, 47] for four fermions 

in a single-j shell, but without an explanation or argument. 

''This does not mean that the 7 = state with nonzero eigenvahie are always lower in energy 
than other / states for Gj = —Sjj'. There are exceptions where other / states appear lower than 
7 = 0. For example, the 7inax state is the lowest for J = 2j — 1. 
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5.4 Schematic interpretation of the empirical approach 

The rationale of the above empirical approach can be seen from the following anal- 
ysis. Although the relation between the eigenvalues and two-body interactions is 
complicated, the eigenvalues are always linear in terms of two-body interactions in 
a "local" space (explained below) of the TBRE. Therefore, instead of studying the 
effects of all two-body matrix elements simultaneously, we can decompose the prob- 
lem into A'^ parts. In each part, we focus on only one interaction matrix element, 
i.e., the entanglements between the two-body matrix elements are neglected. Let 
us take a certain Gj = — 1 and Gj> = {J' ^ J), and diagonalize the Hamiltonian 
of, e.g., fermions in a single-j shell. Suppose that is the lowest eigenvalue. The 
wave function corresponding to Ejp is 

KJ'j 

where (j"~^(i<'7), j^(J')|}j"//9J)'s are the two-body coefficients of fractional parent- 
age (cfp's), which have been widely used in shell model calculations. (3 (7) refers to 
additional quantum numbers needed to define a state of n (or n — 2) fermions with 
total angular momentum / (or K) uniquely. 

Now we introduce a small perturbation by adding e{Gj} to Gj = — 1. Gj = — 1 
and e{G'j} define the "local" space of our TBRE. Let us call it the J subspace. The 
new eigenvalue is approximated in the first order by 

{Ei,)' = Ei, + e'^^^ E [{r-\K^).3V')\mPJ)]"Gr. (44) 

This means that the -E/^ is linear in terms of {Gji} in the J subspace. Because Ej^ 
gives the lowest eigenvalue for the case with Gj = — 1 and others zero, this angular 
momentum / continues to give the lowest eigenvalue in this local space (e can be 
0.2-0.3 in most cases according to our numerical experiments). The full space of the 
TBRE Hamiltonian can be covered mostly by the subspaces defined above. 

Thus the empirical method of Refs. [56, 58] implicitly assumes that a very large 
part of the full space of the TBRE Hamiltonian can be covered by the A^ local 
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subspaces, which are defined by introducing a small perturbation {eGj/} (e is small) 
to a fixed Gj = — 1 (J 7^ J'). For fcrmions in a singlc-j shell, for instance, Gj = — 1 
and {eGj>} define the (J/2 + l)-th local subspace of two-body matrix elements. The 
local subspace with J = is the first subspace, and that with J = 2 the second, 
and so on. 

This philosophy can also be shown from the following numerical experiments: 
Let us take four fermions in a single-j shell with j — ^. In Fig. 12(a) G'j^3^( Jmax = 
16) = —1 and all the other two-body matrix elements Gj> (J' 7^ Jmax) are set to 
be the TBRE, but with a factor e multiphed. One sees that almost all cases of 
the g.s. belong to / = /max — 4j — 6=28 when e is small (say, 0.4). If one uses 
^Jmax(<^max = 16) = 1, then P(/max) ~ 0, which means that the cases of the TBRE 
with Gj^^x < produce almost all the /max g-s. in a single-j shell. Fig. 12(b) shows 
the results for the same system with Gq being —1 and other G/s being the TBRE 
multiplied by e. It is seen that the g.s. is overwhelming for small e. When e = 
and Go — +1, the g.s. probability is also sizable because of the contributions from 
J = 6 and J = 12 (refer to Table VI). 

The above numerical experiments are not trivial. By this procedure one can find 
which interactions, not only monopole pairing, are important in obtaining the g.s. 
dominance. Taking four fermions in a single- j {j — 31/2) shell as an example, the 
g.s. probabihty is ~0.2% if we delete all two-body interactions which produce / = 
g.s. ( J = 0, 6, 8, 12, and 22, refer to the last row of Table VI). This means that the 
g.s. dominance comes essentially from those five interactions. Previously, Johnson et 
al. noticed that the robustness of the g.s. dominance is more or less independent of 
monopole pairing [12, 13]. It was not known, however, whether a certain two-body 
matrix element is essential or partly responsible, and how to find which interactions 
are essential, in producing the g.s. dominance for a given system. 

A shortcoming of the empirical approach of Refs. [56, 58] is as follows: One 
takes one of Gj's is set to be —1 (other Gj's are zero) in each numerical experiment, 
and one finds the angular momentum of the lowest state. The case with Gj = 1 is 
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excluded in the numerical experiments because one would likely obtain degenerate 
levels with the lowest energy for this case. Thus the good consistence between the 
predicted / g.s. probabilities by the empirical approach of Refs. [56, 58] and those 
obtained by diagonalizing the TBRE Hamiltonian indicates that the properties of 
local spaces defined by {Gj — —l + eGj'{J' 7^ J)} more or less represent the features 
of the full space ^. 

5.5 Spin /max g-s. probabilities 

For fermions in a single-j shell, the state with the maximum angular momentum 
(denoted as /max) was found to have a sizable probability to be the g.s. [45, 47]. 
This phenomenon can be explained by the observation that Mi^^^ = 1 always ^. The 
predicted /max g-S- probabilities of fermions in a single-j shell are — j+1/2 ^ 100%) 
independent of all particle numbers (even or odd), which is shown in Fig. 13(a). 

The above argument on /maxg-S. probabilities for single-j shells can be readily 
generalized to many-j shells. Consider, for example, two shells with angular mo- 
menta ji and j2. Following the same logic as was used for a single-j shell, we predict 
that the two angular momenta 1^^^ — /max(i") and /max(i2) have g.s. probabilities 
which are at least as large as x 100%. Here, ImaxUi') (i = 1 or 2) is the largest 
angular momentum of a state constructed from the configuration, and is the 
number of independent two-body interactions in the (ji, j2) shells. In other words, 
we can predict in this way the lower limit for these I^^^g.s. probabilities. 

Figure 13(b) presents I^^^ = /max(jT) and /max(j2) g-S- probabilities. They 
are compared with a simple 1/N plot. Indeed, the predicted lower limit of 1/N for 

^However, some {Gj = l + eGj'{J' ^ J)} local subspaces also contribute to / g.s. probabilities. 
In some cases we find that for an even number of fermions the g.s. probability may not be very 
small even if one deletes all the terms which give / = g.s. with only one of Gj = —1 and other 
Gj- switched off. For example, forn = 4 and j = ^ the P(0)=50.2% if we use the full TBRE 
Hamiltonian, and F(0)=23.1% if we delete Gj's with J = 0,4,8,10 which present / = g.s. if 
these Gj = -1 (refer to Table VI). If we deleted G12 then the P(0) would be 14.2%. Nevertheless, 
numerical experiments provide a very simple recipe to find what interactions are essential for a 
certain / g.s. probability. 

^In Ref. [58], an argument was given why Ej^^^ is the lowest eigenvalue when Gj^^^ = — 1 and 
other Gj's are zero. 
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-^maxS-S- probabilities works very well. It should also be noted that Jg.s. probabilities 
with / very near I'^^^^ are extremely small (less than 1%) in all these examples. 

Figure 13(c) shows the P(/max)'s vs. spin / of bosons, with n ranging from three 
to six. When / is small the agreement is good; when / becomes larger, deviations 
between the values calculated by diagonalizing the TBRE Hamiltonian (denoted as 
P™^^(/niax, 0) those predicted by the above 1/N — 1/(^ + 1) relation appear. 
The -P™^^(/max, are systematically larger than 1/N and increase with n. 

An argument why the behavior of the -P™^^(/max, 0'^ different for bosons and 
fermions, given in Ref. [50], is as follows. As discussed above, the -P(/max) comes 
essentially from a gap produced by the pairing interaction Gl^^^ for bosons with spin 
I, or the pairing interaction Gj^^^ for fermions in a single-j shell, where Lmax = 2i 
and Jmax = 2j — 1, respectively. One can trace the "anomaly" of -P'^^^^(/max) 
back to the gap associated with Gj^^^ or Gl^^^ by using analytical formulas of 
{Ei^^^_2 - -E/^nax)- Here the state with / = /max (-^max - 2) is found to be the 
ground (first excited) state if Gj^^^ or Gl^^^ is -1. We obtain {Ei^^^_2 - Ej^J as 
follows, 

2ln-l 



boson systems : 
fermion systems, n — A : 



Al-1 



105 



+ 



135 



128(4j - 7) 64(4j - 5) 
63 



128(4j-3)' 



r • . r 35 

fermion systems, n = 5 : — - + 



2205 



+ 



5145 



128 2048(4j - 9) 2048(4j - 7) 
1785 189 



2048(4j - 5) 2048(4j - 3) ' 



27 

fermion systems, n — 6 : — — + 

128 



10395 



2835 



8192(4j-ll) 1024(4j-9) 
4725 45 



297 



4096 (4j - 7) 256 (4j - 5) 8192 (4j - 3) ' 



(45) 



One easily sees that the gap for bosons with spin I increases regularly with n at an 
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interval^ 1/2 if / is large, while that for fermions in a single-j shell is much smaller 
(almost one order) in magnitude and comparable for different n and j. For instance, 
the gap is 0.47, 0.39, 0.35 for n =4, 5 and 6 fermions in a j = 15/2 shell, respectively, 
while the gap is 2.03, 2.56, 3.07 for n = 4, 5, 6 bosons with spin / — 7, respectively. 
According to the empirical rule discussed above, a relatively larger gap makes the 
corresponding P{I) larger: the larger the gap is, the larger the corresponding P{I) 
is. It is pointed out again, however, that the the -P™^^(/max,0 follows the 1/N 
prediction when the spins of bosons in the system are small. 

For sd bosons, it was found in Refs. [16, 17] that /max g-S- probabilities are large, 
which can be explained in the same way: Among the two-body matrix elements, 
interactions with C4 = — 1 and others being produce the lowest eigenvalue for 
the /max = 2n state. The predicted I = 2n g.s. probability is 1/A^=1/6=16.7%, 
consistent with that obtained by diagonalizing the TBRE Hamiltonian (~ 15%) 
In sdg-hoson systems, the predicted /max = 4n g.s probabilities is 1/A^ ~ 3.2%, 
where N — 32. The /max g-s. probability that we obtain by diagonalizing the TBRE 
Hamiltonian are 3.3%, 4.2%, 3.3% for n —A, 5, 6, respectively. 

Therefore, for fermions in a single-j shell the 1/N relation works very well for 
both small and large j, and for both small and large n; for fermions in many-j' 
shells, the 1/A^ relation predicts very well the lower hmit for / = I'^^ — /max(jf) 
g.s. probabilities. The P™^^(/max, 0'^ bosons follow the 1/N prediction when 
the spins (/) of the bosons in the system are small. Large deviation from the 1/N — 
l/(/ + 1) relation appears when / is large. 

5.6 Argument based on width 

Two definitions of width were used in the literature. One is defined by gj — 



Di, where Dj is the number of states with angular momentum /. The 



^°Note that the term {s^d^){sd) gives degenerate lowest eigenvalues for many / states when 

Ssdsd is set to be —1 and others arc 0. Therefore, wc use six (instead of seven) as the number of 
independent two-body matrix elements, N. The difference due to this minor modification is very 
small (~ 2%), though. 
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width 's in Table IV are calculated by using this definition. The other is defined 



by CT/ = y {{H — EiY) / Di. In this subsection and Sec. 6.1, the overhead line above 
a certain quantity (such as Ei) means to take the average of this quantity over all 
the states with spin / obtained by one set of the TBRE Hamiltonian, and ( ) outside 
means to take the TBRE average (in our case averaging over 1000 runs). 

As discussed in Sec. 4.1, the argument based on gi is not applicable to cases of 
fermions in a single-j shell. For example, although the / = /max state for fermions 
in a single-j shell has the the largest width gi, P(/max) is smaller than P(0) for all 
j > 7/2. Here -P(/max) was shown to be + \) in Sec. 5.5. The width gi does 
not have a clear relation or correlation with P{I) for fermions in a single- j shell. In 
Ref. [17], Bijker and Prank found a similar result for sd bosons. 

Very recently, Papenbrock and Weidenmueller refined the width [gj) argument 
in Ref. [70]. They derived the distribution of and the correlation between the 
gi''s. They empirically obtained that the maximum of I-E//3I — rjgi (/ fixed), where 
the value r/ is a constant and is called "scaling factor". Por a single-j shell with 
n — Q and j — 19/2, r/ ~ 1.8 for small and medium 7, and ~ 1 for 7 ~ 7inax- 
Interestingly, they found that the 7-'(7) ~ probability for rigi to be maximal for six 
or eight fermions in a j = 19/2 shell. It was speculated in Ref. [70] that similar 
considerations would also apply to other many-body systems. 

In Refs. [71, 72] we tried to understand P(7)'s by combining the property of 
width (J/ and the probability for Ej to be the lowest energy (this probability is 
denoted as 'P{I) and will be be discussed in details in Sec. 6). The idea of Refs. 
[71, 72] is as follows: P(7min) and P(7max) are always large, which can be traced 
back to the quasi-randomness of two-body coefficients of fractional parentage (as 
discussed in Ref. [73]). On the other hand, one should be aware that there is only 
one state with 7 = Ima^i while the number of 7 = 7niin=0 states for an even number 
of fermions in a single-j shell is usually larger than one. Some of these 7 = states 
are pushed down far from their average energy. 
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We define 



a 



J 



n{n — 1) 



T.{r^K^)fimrm{r\K^mj)\}rp'i) , 



2 



and ajp — ajpp. The coefficients {j'^~'^{K^)f{J)\}j'^(5I) are the two-body cfp's. 
The aip such defined is consistent with that used in Sec. 3 and Sec. 4. 

For the case of fermions in a single-j shell, Refs. [71, 72] interpreted large P(/mm) 
and P(/max) in terms of the large fluctuations of aj = Yj/sf^'iis/Di, and large ai for 
1 = states in terms of the statistical point of view [71, 72]. 

One thus expects that the probability of a O"*" state to be the ground state is 
larger than that of the / = /max state. This is a reasonable argument for the 
g.s. dominance for a system with an even number of fermions, although it is at a 
qualitative level. 

Along the fine of Refs. [71, 72], the case of four fermions in a single-j shell can 
be further elucidated by evaluating ct/'s for 7 = 0, 2 and 4. These cr's are found 
empirically to take the largest values. It was shown in Ref. [71] that 



^^One should be aware that this definition of aj is different from gj in Table IV or that in Ref. 
[75], where the energy ccntroid is not taken into account. 

'^^Here one needs to know the number of nonzero eigenvalues Ej^p for states with the angular 

momentum I when Gj = —6jji. For 7 = there is always one nonzero eig envalue E^'^^^ corre- 
sponding to Gj = —Sjj' (refer to Sec. 5.3). For the cases of 7 = 2 and 4 we give the numbers 
of nonzero eigenvalues Ej^p {i is the index for the nonzero eigenvalues) for Gj = —5jji without 
details: 



! = {{Elf) - {{Ej)") 



7 = 2 




1 if J = 
3 if < J < J, 



max 



max 




-2 



max 



-2 
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with 



((^/)') = EE(«w)V(^/) 

J 13,13' 

and 

((^)'> = e(e<) KD:f. 

The above expressions for {{Ej)"^) and ((-E/) ) can be further simphfied in terms of 
non-zero eigenvalues Ej^p (z is the index for the nonzero eigenvalues, the number 
of Ej^f^ is usually much smaller than Dj) for Gji = —Sjjr. 

((^)=eeK,H'/(^^) (46) 

J i 

and 

((^)') = E(E^/?)'/(^/f ■ (47) 

Although it is oversimphfied to assume that all non-zero eigenvalues Ej^f^s are 
equal, it is very instructive to estimate the ct/'s by using this assumption. One can 
then obtain 

^?-(i + ^)EK,H'/^^-36/(i + ^)- 

By using the analytical expressions of the number of states with angular momentum 
/ for four fermions in a single-j shell [48], one finally obtains 

for/ = 0: (7|-12-36/(j + ^) , 
for / = 2 : a? ~ 8 - 36/(j + ^) , 
for / = 4 : (7^ ~ 7 - 36/(j + ^) . 



^^Let us take / = as an example. One has ^■'^ ~ &D^P /{j + i) ~ 2, where d'^P ~ j/3. This 



2 

gives ctq ^ \/l2 in the large j limit. 
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In the large j limit, these ai saturate at ~ \/l2, \/8 and \/7, respectively. For 
j = Y' '^i obtained by 1000 runs of the TBRE Hamiltonian is 3.52, 3.16 and 
2.97 for / = 0, 2 and 4, respectively, which are reasonably consistent with the above 
estimations, ^12 ~ 3.46, ^8 ~ 2.82, and ^/l ~ 2.64. 

One therefore sees that the width cx/ for / = states is larger than those of states 
with other /'s. Because P(0) is large, the P(0) is expected to be even larger due to 
the large gq. This leads to the spin g. s. dominance for the cases of four fermions 
in a single-j shell, and similarly for the case of four bosons with spin I. 

For cases with very large dimensions (larger than 100) and large enough particle 
numbers, Zuker and collaborators presented a formula [74] to obtain the lower bound 
for the energy of states with angular momentum I: Ej — ^\nDi/\n2cri. In Ref. [62] 
Velazquez and Zuker made an effort to relate this lower bound to the problem of 
the g.s. dominance in the presence of random interactions. One of the conclusions 
of Ref. [62] is that the width aj, rather than the energy centroid Ei, plays a crucial 
role for the lower bound energy. However, the origin of large ctq is not yet clear, 
except the case of four fermions in a single-j shell as discussed above. 

5.7 Relation between g.s. wavefunctions of systems with 
mass number differing by two 

In Ref. [13] it was reported that the pairing phenomenon seems to be favored 
simply as a consequence of the two-body nature of the interaction. The "pairing" 
here means that there is a large matrix element of the S pair annihilation operator 
between the ground states of a n fermion system and a n + 2 fermion system in the 
same shell. This seems to suggest that the spin ground states obtained by using 
the TBRE Hamiltonian are, to a large extent, S'-pair condensation. In order to test 
this hypothesis, Johnson et al. followed the example of generalized seniority and 
considered the general pair creation and annihilation operators — J^jCj^p 
S = EjCjSj where S] = (44)^°^ and S = {3^)1 The coefficients is given by 
{n\Sj\n + 2) for samplings with g.s. for both the n and n + 2 fermion systems. (It 
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is noted that the S pair such determined is different from that determined in Sec. 
6.2.2). The pair-transfer fractional collectivity from the g.s. of n + 2 particles to 
that of n particles is defined as follows: 



Thus the 7 = ground states are condensates of S pairs if /transfer = 1- They 
showed that /transfer's between 7 = ground state of n and that of n + 2 are close 
to /transfer = 1 iustcad of 0, Suggesting a correlation of the pairing-like condensates. 
The pair-transfer fractional collectivity may be defined in another form: 



In order to investigate the correlation between the g.s. of a system with n 
fermions and that with n + 2 fermions in the same shell, the case of fermions in a 
single-j shell (where the seniority quantum number v is well defined), was checked 
in Ref. [58]. Below a few examples are discussed. 

The simplest case is four and six fermions in the j = 11/2 shell. The g.s. 
probability for n = 4 and 6 is 41.2% and 66.4%, respectively. Among 1000 sets of 
the TBRE Hamiltonian, 364 sets give g.s. both for n = 4 and 6 simultaneously. 
Namely, the TBRE Hamiltonian by which the ground state has spin 7 = for 
n — A has an extremely large probabihty (around 90%) to produce the 7 = ground 
states also for n — 6. The /transfer's are in most cases around 0.8-0.9 and /transfer's 
are typically around 0.9-1.0. This means that the S annihilation operator takes 
(approximately) the g.s. of six fermions to that of four fermions. 

Now let us take a larger value of j, i.e. j = 15/2, which is good enough for 
our discussion. The g.s. probability for n = 4, 6 and 8 is 50.2%, 68.2% and 
32.1%, respectively. Among 1000 sets of the TBRE Hamiltonian, We found 310 sets 
which produce g.s. simultaneously for n =4, 6 and 8. Considering this 31% and 
P(0) = 32.1% which produce g.s. for n = 8, we can say that almost all those 
TBRE Hamiltonian which produces g.s. for n — S, also produce g.s. for n = 4 



transfer — 



((w,7 = 0g.s.|^|n + 2,7 = 0g.s.))' 
{n, I ^ g.s.\S^ S\n, I ^ g.s.) 



(48) 



transfer 




(49) 
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and 6. It should be noted that the difference of g.s. probabihties among those for 
n = 4, 6, and 8 arc large. As discussed above, the g.s. probability of six fermions 
is 68.2% while that of eight fermions is 32.1%, which means that more than 50% of 
the g.s. for n = 6 are not related to the chain in which the g.s. of n fermions 
can be obtained by annihilating one S pair from that of (n + 2) fermions. 

Figure 14 shows the distribution of seniority in the g.s. which are obtained by 
using the TBRE Hamiltonian, for a few cases of four and six fermions in a single-j 
shell. Low seniority components in the wavefunctions of these g.s. are not favored 
at all. This means that the contribution to the total g.s. beyond the seniority 
chain described in Ref. [13] is more important in the g.s. of these systems. 

Based on the above discussions, we conclude that a chain of angular momentum 
zero ground states, which were suggested in Ref. [13] to be linked (approximately) 
by the S pair operator, can be also seen frequently in systems with even numbers 
of fermions in a single- j shell for small j. However, this chain covers only a part of 
the g.s., because the contribution beyond this chain can be more important, and 
this link becomes weak for large j. 

As for fermions in many-j shells, recent calculations [61] showed that the seniority 
distribution in the g.s. is very complicated. In the sd shell systems low seniority 
states do not dominate in the spin zero ground states. 

5.8 Other results 

One alluring suggestion on the origin of the g.s. dominance is time reversal T 
invariance. Because the time reversal invariance plays a key role in the formation of 
0"*" pairs in the ground states of even-even nuclei, one expects that this invariance 
may imply a built-in favoring of / = ground states in the presence of the TBRE 
Hamiltonian. To see whether or not this is true, Bijker, Prank and Pittel [75] 
analyzed a system of identical nucleons in the sd shell. They took the Gaussian 
unitary ensemble for two-body interactions, for which the time-reversal invariance 
does not hold, rather than the TBRE. The outcome of their calculations was that 
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the g.s. dominance is even more pronounced, which suggests that the time-reversal 
invariancc is not the origin of the g.s. dominance ^'^ Instead, these authors showed 
that for the cases which they checked, / = states have a larger width than other 
/ 7^ states. 

Another view was suggested by Drozdz and Wojcik in Ref. [76]. They found that 
the non-zero off-diagonal matrix elements of all GjS for / = states have a wider 
distribution on average than those of / 7^ states, and thus the 7 = states are 
expected to spread over a broader energy interval even though the number of states 
is usually much smaller than for the / 7^ cases. However, one can not conclude 
that the g.s. dominance arises just from this phenomenon. As pointed out by 
the same authors [76], the g.s. dominance results from an interplay between the 
diagonal and off-diagonal matrix elements. 

In Ref. [77], Kaplan and Papcnbrock studied the structure of cigcnstatcs for 
many-body fermion systems in the presence of the TBRE Hamiltonian. They found 
that near the edge of the spectrum, wave function intensities of the TBRE Hamil- 
tonian exhibit fluctuations which deviate significantly from the expectations of the 
random matrix theory. A simple formula was given which relates these fluctuations 
to the fluctuations for the TBRE Hamiltonian. The possible connection between 
this deviation and the g.s. dominance in fermion systems with even numbers of 
particles is unclear. 

Because large energy gaps were found in Refs. [12, 13], between the g.s. and 

excited levels, one can ask whether there exist certain universal features of these 

gaps. In Ref. [78], Santos, Kusnczov and Jacquod defined the energy gap as follows: 

For sp {sd) boson systems, it is defined by the energy difference between the first 

1~ (2+) state and the / = ground state energy; For an even number of fermions 

^''in Ref. [62] Velazquez and Zuker claimed that the time reversal invariance suggests the g. 
s. dominance, though it does not imply it. Namely, the g.s. dominance must be associated 
to some general cause, and the time reversal invariance is a good candidate, whose influence can 
be detected through the abundance of self-conjugate T pairs. For fermions in a single-j shell, 
considerably large components of g.s. are given by low seniority pairs, in particular for cases 
with j not large. This alluring suggestion deserves further studies. 
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in a single-j shell, it is defined by the energy difference between the first l"*" state 

and the / = ground state (after a re-scaling). They showed that the distribution 
of the gap such defined is robust and may be helpful in understanding the g.s. 
dominance. 

Spin-1/2 fermions (i.e., without orbital angular momentum) in the presence of 
random interactions were studied recently in Refs. [79, 80, 81]. Jacquod and Stone 
[79] derived a formula for fixed-/ (total spin of the system) widths for the TBRE 
Hamiltonian by using a heuristic argument, while Kota and Kar [80] used a group 
theoretical approach coupled with a so-called binary correlation approximation. Ka- 
plan et al. [81] showed that one has to go beyond the width and consider a so-called 
excess parameter which gives deviations from Gaussian distributions. Adding this 
correction, one sees the smallest / to be lowest in energy for spin-1/2 fermions. 

5.9 Summary of this Section 

In this Section we have reviewed the results for complicated systems. We first 
explained and applied an empirical rule to predict P(/)'s. This rule was found 
to work for fermions in a single-j shell or in many-j shells, with even numbers of 
particles or odd numbers of particles. The same rule works for bosons. 

The empirical rule also means that the g.s. dominance is essentially given by 
two-body matrix elements which produce 1=0 g.s. when one of two-body matrix 
elements is —1 and all others are zero. For fermions in a single-j shell or bosons 
with spin /, the origin of the /max g-S- probability is clarified: The large P(/max) 
is essentially contributed by an attractive Gj^^^. The P(/max) would be close to 
zero if one deleted Gj^^^ for fermions or Gl^^^ for bosons. The simple relation 
-P(-^max) ~ 1/-^ works vcry well for fermions and also for bosons with small I. 

The disadvantage of this approach is that one must diagonalize the Hamiltonian 
under the requirement that one of two-body matrix elements is —1. Therefore, this 
interpretation is not so transparent. 

We further reviewed our results presented in Refs. [50, 68], where an argument 
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for the g.s. dominance and the regular staggering of P(0) versus j was given for 
four fermions in a single-j shell and four bosons with spin /, based on this empirical 
rule. This simple argument was found to be restricted to n = 4 with single-j or 
single-/ shell and to a large j (/) shell with n = 6. It is difficult to "generalize" 
to other cases. The disadvantage of this argument is that it does not provide a 
quantitative evaluation of P(/)'s. 

There were a few interesting efforts to relate [62] the g.s. dominance to the time 
reversal invariance of the Hamiltonian. None of these alluring arguments provides 
us a good description of / g.s. probability. In contrast, it was found in Ref. [75] 
that a Hamiltonian which breaks the time reversal symmetry enhances the P{0). 

There were also arguments [71, 72] of the g.s. dominance in terms of the large 
value of width (Tq defined by \J {{H — £'/=o)^)/-D/=o- First, it was realized [71, 72, 73] 
that the probability 'P(O) for energy centroids with / ~ and / ~ /max to be the 
lowest is large. Thus -P(O) is large. Second, it can be shown schematically, at least 
for the case of four fermions in a single- j' shell or four bosons with spin that the 
7 = states have the largest width (Jo, based on which one expects that the P(0) 
should be even larger. This argument is interesting but it is unknown yet how to 
relate P{I) to V{I) and aj in an explicit way. 

We also reviewed a regularity of the so-called pair-transfer fractional collectivity 
/transfer, which is defined by matrix elements (after normalization) of the S pair 
operator between the g.s. of systems with fermion numbers differing by two. It 
was found in Ref. [13] that the /transfer values are large (larger than 0.5 in most 
cases) for even numbers of valence nucleons in the sd shell. In Ref. [58] it was found 
that for the case of single-j shell with small j the /transfer values are almost 1.0. One 
should be aware, on the other hand, that the chain of the g.s. for systems with 
mass number differing by two covers only one part of the g.s., and /transfer is not 
large when j is large. 

The results in Refs. [76, 77, 78] might be also helpful towards understanding of 
the g.s. dominance. Interesting results include an observation for a robust energy 
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gap between the g.s. and some specific excited states [78], an observation of a 
wider distribution of off-diagonal matrix elements for / = states [76], discussion of 
structure of eigenstates obtained by random interactions [77]. However, the resuhs 
of Refs. [76, 77, 78] are far from apphcations to detailed analysis of P(/)'s. For 
instance, these authors were unable to give any predictions of P(0)'s for the systems 
that they studied. The analysis in Refs. [79, 80, 81] are restricted to systems with 
spin- 1/2 systems. 
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6 Average energies, collectivity and yrast state 
spin 



In previous Sections we discussed the regularities of the ground states in the presence 
of the TBRE Hamiltonian. In particular we focused on P(/)'s and the approaches 
to predict them. 

In this Section we shall the behavior of average energies, collective motion and 
normal ordering of yrast state spin of many-body systems in the presence of the 
TBRE Hamiltonian. We shall also discuss some results on one of the original ques- 
tions, how arbitrary an interaction can be in order to reproduce the global regular- 
ities of atomic nuclei. 

6.1 Behavior of average energies 

There were a few studies addressing average energies. It should be noted that the 
meaning of "average energies" used by different authors can be different, and one 
should be careful about their definition. 

6.1.1 Definitions of average energies 

Let us begin with average energies defined in Refs. [12, 82], where the authors 
averaged energies of the yrast states with even values of /. Ref. [12] showed an 
indication of /(/ -|- 1) behavior for these average energies, which was called "non- 
collective" rotation in nuclear spectroscopy. The authors of Ref. [82] investigated 
the transition from the realistic two-body interaction to purely random ones (also 
refer to Sec. 6.3), and showed that the relative ordering among the energies of the 
yrast states averaged over the TBRE survives with exceptions of low /, but their 
relative separations change significantly when two-body interactions change from 
realistic to random ones. 

In Refs. [45, 46] Mulhall et al. chose a subset which gives 7 = g.s. or 7 = /max 
g.s. among 1000 runs of the TBRE Hamiltonian. Then they averaged the energies 
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of the yrast states for this subset of the ensemble, and found that these average 
energies of the yrast states such defined follow a parabolic function of / with a 
strong odd-even effect. 

Let us change here the definition of average energies. From now on, we take 
the centroids (denoted as Ej) of all eigenvalues of states with spin I as the average 
energies, as defined in Sec. 5.6. One sees that Ej is a linear combination of Gjs: 

E'i = YyiGj, (50) 
J 

where aj is obtained by averaging a/^ over all /9's. One can apply the empirical 
approach of Sec. 4.2 to predict the probability (denoted as V{I)) for Ej to be the 
lowest energy. 

The motivation to investigate Ei is that its behavior is much simpler than that 
of the Ejps, and we expect a more transparent explanation of its behavior. This 
explanation might be very helpful to understand the original problems such as the 
g.s. dominance of many-body systems in the presence of the TBRE Hamiltonian, 
as we have discussed in Sec. 5.6. 

6.1.2 The probability for Ej to be the lowest energy 

One should be aware that a displacement of the TBRE produces only a constant 
shift to all E^/'s. This can be easily seen from Eq. (50). A displacement TBRE 
with c defined by Eq. (17) gives the same Ei plus a shift \n{n — l)c, because 
Y.joij — \n{n — 1). Thus all regularities of -B/'s, including probabilities V{iys, are 
robust regardless of the displacement, while P{I) is very sensitive to the displacement 
(refer to Sec. 3.5). 

In Fig. 15 we plot 'P(/)'s for a few different systems: four fermions in a single- j 
shell {j = y), six fermions in two-j (2ji, 2j2) = (ll,7) shells, six sdg bosons, and five 
fermions in a single-j shell (j = |). They are typical examples among the many 
cases that we have checked: four, five, and six valence fermions in a single- j' shell 
up to J = Y, both even and odd numbers of fermions (n = 4 to 9) in two-j shells 
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with (2ji,2j2) = (7,5), (11,3), (11, 5), (11,7), (11,9), (13,9), c^-boson systems with na 
changing from 3 to 45, srf-boson systems with n changing from 4 to 17, and sdg- 
boson systems with n = 4, 5, and 6. One sees that the P(/)'s are large if / ~ /mm 
or /max, and are close to zero otherwise. This feature holds for all the cases that we 
have checked. 

Now let us describe the 7^(/)'s by using the empirical approach in Sec. 4.2. As a 
specific example, we discuss four fermions in j = | shell. We predict the ■P(/)'s by 
using integrals similar to Eq. (20), without diagonalizing the TBRE Hamiltonian. 
The predicted ■P(/)'s for four fermions in j = | shell are listed in the column "predl." 
in Table VIII. One sees that the V{I) is large if Ej has one or more aj which are 
the largest (or the smallest) for different /'s. The P(/)'s predicted by Eq. (24) are 
listed in the column "pred2." of Table VIII. In Eq. (24) A/} is the number of times 
for a/ to be either the smallest or the largest for each /, and N„i = 2N —1 = 9. 
The P(/)'s obtained by using 1000 sets of the TBRE Hamiltonian are given in the 
column "TBRE". It is seen that the two predicted P(/)'s are reasonably consistent 
with those obtained by using the TBRE Hamiltonian. 

The origin of the fact that ■P(/)'s are large only if / ~ /min or / ~ /max was 
argued in Refs. [45, 73] based on the assumption that the two-particle cfp's are 
randomly distributed, and in Ref. [80] based on group symmetries of the TBRE. 
We shall discuss the formula of Ej derived in Refs. [45, 80] in Sec. 6.1.4. 

6.1.3 The /(/ + 1) behavior of 'E~i 

Examining the ordering of the average energies Ej^s with respect to /, one sees 
that, when the spin / of the lowest Ej is ~ /min (-^max) OHC sct of the TBRE 
Hamiltonian, Ej given by the same set of interactions increases (decreases) with / 
in most cases, and is proportional to /(/ + 1) on average. One should be aware 
that both even and odd values of / are included here. Although the authors of 
Refs. [12, 82] discussed the /(/+ 1) behavior of average energies, their definition of 
averaging differs from our Ej, and they were interested only in even / values (see 
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Sec. 6.1.1 for their definition). Also, one should not confuse the Ej^s with Bethe's 
expression for the level densities [27] which is is based on the Fermi gas approach 

Let (-E'/)min {{Ei)ma,x) bc a quantity obtained by averaging the energies Ej over 
the cases of Ei^^i^.^ (-^/~/max) being the lowest. We find that both (£^/)min and 
{Ei)raax are proportional to /(J + 1), similar to "rotational" spectra. 

In Fig. 16 we show (£'7)min vs. /(/ + 1) for twenty d bosons, ten sd bosons, 
four fermions in a single-j {j—17/2) shell, and a system with four fermions in two-j 
shells. For the sake of simplicity, we introduce the "moment of inertia" J, defined 
by the optimal coefficient as (-E'/)min = 2^-^(-^ + l)- Below we do not discuss (-E'7)max 
but note that {Ej)^^ ~ ^ [/max(-^max + 1) - + 1)], where the J for (Ej)^^ 
and that for (-E'7)min are approximately equal for all cases that we have checked. 

An empirical relation between J' and j is summarized in Fig. 17. For fermions 
in a single-j shell, d boson systems and sd boson systems, J' is fitted by a trajectory 
of the form ~ l-42j (we take j = 2 for d boson and sd boson systems); for sdg 
bosons and fermions in two-j shells, there seems a slight shift from the trajectory 
y/j ~ 1.42j, where f = + jl 

Another feature of V{iys is that the 7'(/max)'s are always quite "stable" (about 
28-35%), while the V{Iinax — 2)'s of fermions in a single-j shell and boson systems, 
and the V{Imax. — l)'s of fermions in many-j shells, are drastically smaller than the 
'^(-^max)'s, though still sizablc (~ 5 — 15%). Moreover, there may be 2 or 3 sizable 
P(/)'s for the cases of / ~ /min, and the P(/min)'s are not always larger than the 
other P(J)'s (with / /mm)- 

In Ref. [73], an argument of the above asymmetry for V{I) was given in terms of 
the fiuctuations of Ej: The Ej^s are proportional to /(/ -|- 1) with some deviations 
in each set of the TBRE Hamiltonian. Because (-Bz+i — eA is small if / is small 



^^The 7(7 + 1) behavior of Ej in Bethe's expression for the level densities and that of {Ej)^in (or 
(7^/)max) discussed in this paper are completely different. For example, J in Bethe's expression 
changes with particle number, but in the present context it is not sensitive to the particle number 
n but to the orbits of the shell, as shown in Ref. [73]. Furthermore, systems which show the 
7(7 + 1) behavior of (7?/)inin can be very simple, and those described by Bethe's expression require 
complexity in energy levels so that one needs a statistical approach. 
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and large if / is large, the possibilities to change the order of the Ei^s due to the 
deviations when / ~ /min arc much larger than when / /max- This explains why 
there are 2 or 3 sizable and comparable P(/)'s with / ~ Jmin but only one large 

with / = J^ax, P(/max " 1) < V{Ira^), and P(/,„ax " 2) < P(/max). 

The relation (-E/)min ~ -^(-^+ 1) was also discussed in Rcf. [73] by assuming that 
the two-body cfp's are randomly distributed. The analysis in Ref. [73] showed that 
such an assumption can explain the /(/ + 1) behavior of (£^/)min- 

The J for fermions in a single-j shell can be evaluated in the following way. 
Let us assume that {Ei=i^^)^ir, ~ Q;/^^(G'j^,Jmin by neglecting contributions to 
Ej^^^ = Ej^^^ from other Gj's. If one uses {Gj^^J — 0.7, as found empirically in 
Refs. [46, 73] for four fermions in all single-j shells, one obtains 

<rAGj^J^n ~ 0.7a/- ~ ^/ma.(/ma. + l) = ^ (4j - 6) (4j - 5) ~ ^f. (51) 

It was shown in Ref. [58] that the a/™^ of n = 4 saturates quickly at y when j 
increases. Then the left hand side ~ y x 0.7 ~ 2.54. One finds that \^ ~ 1.77j, 
which is shghtly larger than the \^ obtained in Fig. 17 (where \^ ~ l-42j). The 
J' in Eq. (51) is over-estimated, because the contributions from and Gj^^^-i 

were neglected. 

6.1.4 The formulas by Mulhall et al. and Kota et al. 

For fermions in a single-j shell, Mulhall, Volya and Zelevinsky derived a formula 

in Refs. [45, 46] for the expectation value of energy of spin / states by minimizing 

the ground state energy and assuming the statistical point of view for the angular 

momentum couplings in an n-body system ("geometric chaoticity" called by these 

authors) . The energy such obtained is actually the centroid of eigenvalues of spin I 

states namely, the average energy Ei defined by Eq. (50). Their final result was 

^^The energy such obtained was interpreted in Refs. [45, 46] to be the lowest spin / energy. 
However, it was later pointed out in Refs. [58, 80, 70, 73] that it should be -B/, the energy centroid 
of the -E//3's. 
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written as follows, 

+ /(/ + 1) S(2J + i) ^;^i^^l, Gj + 0(P{I + If). (52) 

where 0{P{I + lY) refers to higher I terms which seem to be negligible. The first 
term of this formula is a constant which is independent of /. The second term is 
proportional to /(/ + 1). However, this term does not explain ^/J ~ l-42j (refer 
to Fig. 17). Let us take J ~ 2j and Gj ~ 0.7 for the coefficient of this term, 
then one sees that the second term ~ + !)• This is different from our result: 

~ Ej~ + (wc used J ~ 2j2 here). 

In Ref. [80], Kota and Kar obtained Eq. (52) for Ei by resorting to the group 
structure of U{2j + 1) D 0(3) for n fermions in a single-j shell. They also pointed 
out that the use of the cranking approximation and the Fermi-Dirac occupancies of 
particles in Refs. [45, 46] are equivalent to the approach by using group symmetries 
of the TBRE as in Ref. [80]. 

Although the approach of Refs. [45, 46] does not produce the same results 
as those in Refs. [71, 72, 73], the concept of geometric chaoticity introduced in 
Refs. [45, 46] is essentially the same as the randomness of two-body coefficients of 
fractional parentage introduced in Refs. [71, 72, 73]. Both of these two concepts are 
based on the complexity of angular momentum couplings of n particles. 

6.1.5 A short summary 

As a short summary of this subsection, let us first repeat two robust regularities of 
many-body systems interacting by the TBRE Hamiltonian: 

1. The V{iys are large if and only if 7 ~ /min or /max- 

2. The / dependence of the (-E'/)mm and (-E'/)max is roughly /(/ + 1). 

These regularities have been argued in terms of the statistical distribution (geo- 
metric chaoticity) of the two-particle cfp's. 
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We emphasize here that the P(/)'s, which are related to the randomness of 
two-particle cfp's or another appellation of "geometric chaoticity" for the angular 
momentum couplings, and the P(/)'s discussed in the last Section, are different 
quantities. For examples, P{j) is always sizable for an odd number of particles, but 
V{j) is close to zero in this case; P(/max) is not large (~ l/(j + |)) for fermions in 
a large single-j shell, but 7^(/max) ~ 38%. In particular, P(/)'s are independent of 
the displacement c in the displaced TBRE, but -P(/)'s are sensitive to c. We should 
not mix up these two probabilities. 

This shows that / g.s. probabihties (and g.s. dominance) cannot be explained 
by the geometric chaoticity. The role of variances or the width of states with spin I 
is important and complicated. This was stressed in Refs. [71, 72, 73, 80]. 

6.2 Collective motion in the presence of random interac- 
tions 

The structure of levels of many-body systems in the presence of the TBRE Hamil- 
tonian is another interesting topic. In Refs. [16, 17] it was shown that both vibra- 
tional and rotational features arise in the IBM [32] in the presence of the TBRE 
Hamiltonian. In contrast, as pointed out by many authors, rotational behavior does 
not generically arise in fermion systems when their interactions are random. It 
was suggested, therefore, that some constraints should be imposed on the random 
Hamiltonian to obtain a generic rotational behavior in fermion systems [85] . 

6.2.1 Vibration and rotation in the vibron model and the IBM 

The occurrence of both vibrational and rotational band structure within the frame- 
works of the sd IBM and the vibron model was elegantly demonstrated in Refs. 
[16, 17]. 

The following ratio 

R = - E,^)/{E,^ - E,^) (53) 
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has been used as collective indicator, and extensively investigated for medium-heavy 
even-even nuclei [86]. R has the characteristic value of 2 for vibrational systems, 
and 10/3 for rotational systems. 

For the vibron model [37], the collective indicator R is given by 

^=(^2+-^o+)/(^i+-^o+) , (54) 

which is equal to 2 in the vibrational limit and 3 in the rotational limit. Because 
the results for sp bosons and those for sd bosons are similar, we give below only the 
results for so? bosons. 

The R values were calculated by using the shell model by Johnson [85]. He 
found that the i? ~ 1.5 with a broad distribution for ^^O ranging from to ~ 3 
by diagonalizing the TBRE Hamiltonian. This means that the rotational band 
structure does not appear for fermions interacting by two-body random forces. 

In Ref. [16] Bijker and Prank calculated low-lying states within the sd IBM by 
using the TBRE Hamiltonian. They obtained that about 60% among 1000 sets of 
the TBRE produce / = ground states. The same sets, which produce 0^ ground 
states, were used to calculate i?2+ and E^+. These and as well as Eq+ 
determine R values in Eq. (53). The i?'s such obtained are larger than 1 for most 
of them. 

Figure 18 plots the distribution of R with boson number n—?>, 6, 10, and 16 
(taken from Ref. [16]). One sees that two sharp peaks emerge at it! ~ 2.0 and 3.3 
when the boson number n increases from 3 to 10. 

The E2 transition rates also take characteristic values for vibrational and rota- 
tional modes. For example, ^^^'\~^\\ take the value of in the vibrational limit 

^ ' B(E2,2^^0j ) n-l 

and the value of in the rotational limit of the IBM [29, 30, 31, 32]. When 

boson number n is large, these characteristic values in the IBM are equal to those 
in the Bohr-Mottelson model: 2 for vibration and I (2o'2ojoo) 1^ ~ Y ^'-'^ rotation. 

A correlation between R and ^^^^'^+^^+1 for n = 16 is shown in Fig. 19 (taken 
from Ref. [16]), where one-body parameter in Eq. (11) was also included as a 
random parameter together with the TBRE. 
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The remarkable rotational peak obtained in Refs. [16, 17] suggested that the key 
for obtaining a rotational peak from a random shell model Hamiltonian might be 
to restrict the space to a collective subspace built from the lowest S and D pairs, 
because these SD pairs are the objects which are approximated by s and d bosons 
in the IBM [32]. To see whether this is true, calculations in the truncated SD-pair 
space were carried out in Ref. [87] by using a general TBRE Hamiltonian defined 
by Eq. (2). 

The calculations of Ref. [87] were performed in the following procedures: One 
first selects only random interactions which produce O''" ground states and 2+ first 
excited states in two-particle systems. Then it is assumed that the collective S and 
D pairs are given by these 0"^ and 2+ wavefunctions of the two-body systems. One 
proceeds to calculate the spectra for six-particle systems by using the same random 
two-body interactions. 

Figure 20 shows the distribution of R values in the S'D-pair truncated shell model 
space for six identical nucleons in the sd, pf and sdg fermion shells, respectively. 
One sees that the distribution of the R values in the sd shell within an SD-pau 
subspace is similar to that obtained in the full shell model space [85] - a broad 
distribution extending to it! ~ 1.3. When one goes to larger shells, the distributions 
become sharper, and shift to the right from the sd shell {R ~ 1.3) to the sdg shell 
{R ~ 1.91). Nevertheless, no sharp peak at ~ 3.33 appears. Prom this it is 
concluded that rotational motion does not seem to emerge from SD truncated shell 
model calculations. Statistically, the full shell model space and the SD truncated 
truncated shell model space defined here give essentially the same results for a 
general two-body interaction. 

6.2.2 Rotational spectra in the 5'Z^-pair subspace 

Because rotational motion does not arise from the SD-pair truncated shell model 
space with the TBRE, let us investigate whether it might appear when we use a more 
restrictive Hamiltonian defined by Eq. (5) and take the SD-pair truncation as well. 
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Here we fix Hi (single-particle term) to be zero, and keep only monopolc pairing, 
quadrupole pairing and quadrupolc-quadrupolc interaction, with tiicir strength pa- 
rameters (i.e., Go,G2 and k) generated randomly. Each of them follows a Gaussian 
distribution with an average being zero and a width being one. 

The calculations were carried out in Ref. [87] in the same procedure as the SD- 
pair truncated shell model calculations by using the Hamiltonian of Eq. (2), except 
that the G'j(jij2i3i4)'s are replaced by the strength parameters Go, G2 and k,. Fig. 
21 shows the distribution of R thus calculated for six identical nucleons in the sd 
shell. When all three strengths are treated on the same footing, one arrives at the 
distribution of R shown in Fig. 21(a). In this case, no sharp rotational peak is 
observed. Instead, a peak appears around R ~ 1.3, with a long tail extending to 
i? ~ 3.1. If the J2m QmQm strength parameter k, is artificially enhanced by a factor 
e (> 1), one arrives at the results shown in Figs. 21b-d. As e is increased, i.e., as the 
quadrupole-quadrupole strength is enhanced, a peak at i? ~ 3.1 gradually appears. 
On the other hand, the probability of > 3.1 remains very small. 

As the size of the shell is progressively increased, the peak at ~ 1.3 gradually 
disappears and another peak at ~ 3.3 emerges. This is illustrated in Fig. 22 for 
six identical nucleons in the pf, sdg, pfh, and sdgi shells even with e = 1.0. For a 
large shell, the peak at ~ 3.3 becomes very well pronounced. 

As fingerprints of occurrence of rotational motion, ratios of the E2 transition 
rates of six nucleons in the sdgi shell (cf. Fig. 22 (d) ) are examined. According to 
the Elliott Model [33], the ratio ^(£22+^0+) ^'"^^ ^'^^ ^^^^ shell (this ratio is 1.16 

in the IBM for three sd bosons), and I'-^.'^Zlf) = bSeS^^o}! = bIeS^^o^S = 0- 

It is seen in Fig. 23 that there is a strong correlation between R and the ratios 
of these E2 transition rates. One sees a concentration of points near the coordi- 

--tes (i^, ggg:gi )=(3.3,o), (i^, gg:g:gi )=(3.3,o), (i^, gg:g:g; )=(3.3,o), 

i§ij65T)) = (3-3, 1-3), respectively. 
Based on the above results it was concluded in Ref. [87] that the rotational mo- 
tion is related closely to the form of two-body interactions. In particular, for systems 
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of identical nucleons there must be a strong quadrupolc-quadrupolc component in 
interactions for tlic rotational motion to occur. 

6.2.3 Rotations based on displaced random interactions 

The philosophy that the occurrence of rotational band structure requires the Hamil- 
tonian to have some specific features was also discussed by Velazquez and Zuker 
in Refs. [62, 82]. Since the B(E2) enhancement is not produced by the TBRE 
Hamiltonian, they took a displaced TBRE which is centered at a negative constant 
c, namely, the displaced TBRE of Velazquez and Zuker is attractive on the average. 
This idea can be traced back [88] to the fact that the reahstic interaction in the pf 
shell, such as the KB3 [83], are mostly attractive. In Ref. [62], Velazquez and Zuker 
found a gradual buildup of B(E2) values as well as a gradual increase of -P(O) in the 
(/7/2P3/2) space as c decreases from to —3. R = 3.3 was found to be also clearly 
favored. 

One criticism [64] was that the magnitude of displacement used in Ref. [62] is 
too large. The width of the KB3 matrix elements is close to their average value in 
magnitude. However, the value of |c| used in Ref. [62] is about five times larger 
than the average value of the KB3 matrix elements. When the magnitude of c is 
so large, the B(E2) values are dominated by the large value of c, according to Ref. 
[64]. 

Velazquez and Zuker suggested in Ref. [62] another possibility to obtain rota- 
tional band structure for fermions. Their suggestion is based on the observation 
that the nucleus ^°Ne (with four valence neutrons) has a rotational spectrum while 
the nucleus ^^Ar (with four neutron holes and two proton holes) has a vibrational 
spectrum although the same realistic two-body matrix elements are used for these 
two nuclei in the sd shell. This difference comes possibly from the changes of the 
mean field. They suggested that one would obtain both vibration and rotation in 
one nucleus by randomizing the single-particle energies with the two-body matrix 
elements fixed. 
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6.2.4 A short summary 

To summarize this subsection, we first reviewed the discovery of generic vibration 
and rotation within the sd IBM. Similar results can be seen for the sp bosons. It 
was found that in a truncated SD-pair subspace vibrations arise for the general 
TBRE Hamiltonian but rotations do not. With a restricted Hamiltonian defined 
by Eq. (5) collective rotations appear. Not surprisingly, the quadrupole-quadrupole 
interaction seems to play a key role in obtaining a peak at i? ~ 3.33. 

According to Refs. [62, 82], a negatively displaced TBRE is also able to produce 
a rotational band structure as well as an enhanced B(E2) transition rate. One 
criticism to this suggestion was that the displacement used in Refs. [62, 82] is too 
large. 

It is interesting to discuss why the IBM with the TBRE Hamiltonian is able 
to give rise to rotations, while the shell model truncated to SD pairs cannot. As 
generally believed, the sd IBM is a consequence of quadrupole and pairing correla- 
tions contained in the realistic interaction. Then random interactions among s and 
d bosons already absorb some parts of quadrupole correlations among the realistic 
interaction. There is no inconsistency, therefore, between the results of Refs. [16, 17] 
by using the TBRE within the IBM, those of Ref. [85] by using the TBRE within 
the shell model, and those of Ref. [87] by using both the TBRE and the restricted 
Hamiltonian defined in Eq. (5) within the SD nucleon pair approximation [41, 42]. 

6.3 Normal ordering of spin in the yrast band 

The first study of normal ordering of spin I in the presence of the TBRE Hamiltonian 
was done by Cortes, Haq, and Zuker [89] more than twenty years ago. These authors 
took the case of ^°Ne nucleus and used a Hamiltonian consisting of the Elliott SU(3) 
component [33] and the TBRE component. Namely, the Hamiltonian of Ref. [89] 
was defined by 

H={1- 6)i^Elliott + bHTBRE, (55) 
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where i^EiUott is the EUiott Q-Q force and //tbre refers to the TBRE Hamiltonian. 
The parameter h runs from to 1 in order to study the effect on the normal ordering 
in the spectrum from the noise, the i/xBRE part. 

It should be noted interestingly that the calculations of Ref. [89] almost observed 
the g.s. dominance but the large variances therein prevented the authors from 
drawing this conclusion. The Fig. 2 (8) in Ref. [89] already suggested that the 
yrast spin zero state is lower than those for spin two, four, etc. on average, even 
when h —1 (the pure TBRE Hamiltonian). 

In Ref. [82] Velazquez and collaborators apphed the similar idea of Ref. [89] to 
^^Mg, ^^Ti and ^^Cr. The motivation was to investigate the transition from reahstic 
two-body interaction to purely random ones. The realistic interaction which they 
used is the Kuo-Brown (KB3) interaction [83] for the fp shell and the Wildenthal 
interaction [84] for the sd shell. Random interactions were taken as the TBRE. 
It was found in Ref. [82] that the average energies of yrast states with different 
angular momenta / keep the ordering of the band when the Hamiltonian changes 
from the realistic interaction to the TBRE interactions when h is small, and that the 
probability that the yrast states keep the ordering is quite large (~ 30%) even for 
purely random interactions, suggesting a strong correlation between these states. 

Another relevant result was given in Ref. [12], where it was shown that there is 
an indication of a so-called "non-collective" rotation in the spectrum of the nucleus 
^^Ca, i.e., /(/ + 1) behavior of the yrast energies obtained by averaging over the 
g.s. subset of the TBRE, where / is even. 

6.4 Constraints on random interactions in nuclei 

As discussed earlier, within both the shell model and the IBM, the TBRE Hamilto- 
nian reproduces some features of nuclear properties. However, the realistic interac- 
tions between nucleons are of course not random. This leads to essential differences 
between calculated results by using the realistic interaction and those by using the 
TBRE Hamiltonian. For example, the observed P(0) of even-even nuclei is 100% 
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without exceptions, while for the TBRE Hamiltonian it is typically 30 — 70%. There- 
fore, although some characteristic properties survive as the interaction changes from 
realistic to random, it would be interesting to study, as was asked in Ref. [85], how 
arbitrary a set of interactions can be in order to reproduce realistic nuclear proper- 
ties. In this subsection we shall review results in this context. 

The IBM is a very proper tool to investigate this context due to its simplicity. 
Although the IBM with the TBRE can explain some typical features of quadrupole 
type collectivity, there are certain quantities which cannot be simply explained by 
this approach. In Ref. [90], Kusnezov et al. investigated this question within the 
sd IBM. One of the features discerned in Ref. [90] is that for all nuclei throughout 
the nuclear chart with 6 ~ 18 valence nucleons outside the doubly closed shell the 
experimental values of R show a peak at 2.3 while the IBM with random interactions 
tends to give a peak at ~ 2. This indicates that we need to restrict our random 
interactions. Instead of using the general Hamiltonian for the sd bosons in Eq. (11), 
these authors used a schematic sd boson Hamiltonian, H — e^rid — k,J2m QmQm, 
which is called the extended consistent Q formahsm proposed in Ref. [91]. Here 
Qm = H^s^dM + d^M^ + xEr„im2(-)"^'(2mi2m2|2M)4^d_„,. They obtained 
statistical ranges oi e^/ k, and x values which give R in the range of 2.2-2.4 for seven 
sd bosons. 

A similar study was performed by Zhang and collaborators in Ref. [92] within 
the geometric collective model (GCM) [93]. Their calculations demonstrated that 
the experimental values of R constrain some ratios between the parameters in the 
GCM potential. 

This problem was also studied preliminarily in the shell model. In Ref. [94], 
Horoi et al. studied a system of four protons and four neutrons in the sd shell 
(corresponding to the ^^Mg nucleus) by using random two-body interactions which 
distribute uniformly between —1 and 1. They found that the overlap between the 
I — T — ground state wavefunctions obtained by random interactions and those 
obtained by the reahstic effective interaction is about 0.02 on the average. The B(E2, 
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25*" — > 0]^) values obtained by using random interactions are typically one order of 
magnitude smaller than those obtained by the realistic interaction, indicating that 
collectivity of those states calculated by using random two-body random interactions 
is not as strong as that obtained by using the realistic two-body interaction. 

More extensive comparison between the collectivity arising from the TBRE shell 
model Hamiltonian and experimental data is necessary but difficult because of the 
huge dimensionality of the shell model. Due to this difficulty, little has been known 
about the constraints on random interactions of the shell model Hamiltonian so far. 

6.5 Summary of this Section 

In this Section we first discussed the behavior of average energies, and the occur- 
rence of vibrational and rotational band structure, for many-body systems in the 
presence of random two-body interactions. Then we discussed normal ordering of 
yrast spins calculated by using the TBRE Hamiltonian. We also reviewed the re- 
sults of constraints on the random Hamiltonian in order to reproduce the global 
properties (such as the distribution of R) of realistic nuclei. 

Energy centroids Ej were discussed in many papers [45, 46, 62, 80, 71, 72]. 
It was found that the probabihty for Ej to be the lowest is large if 7 ~ /min or 
~ -^max- We thus divide the TBRE into two subsets, one of which gives Ej^i^.^ the 
lowest energy, and the other of which gives Ej^j^^^^^ the lowest energy. The {Ej)^[n 
{{Ei)mm), obtained by averaging the Ej over the Ei^j^.^ (-^/~/max) subset, the 
behavior. These features can be explained in terms of the quasi-randomness of two- 
body coefficients of fractional parentage. 

The occurrence of vibrational and rotational structure for the sp- and sd-boson 
systems with the TBRE was discovered and discussed in Refs. [16, 17]. However, 
the rotational motion does not arise in fermion systems if one takes the general 
TBRE Hamiltonian. Additional requirements are thus necessary to obtain a ro- 
tational band structure for fermion systems. Two kinds of random Hamiltonians 
have been discussed so far to obtain rotational band structure for fermion systems: 
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a Hamiltonian by taking a displaced TBRE with an attractive average [62], and a 
restricted Hamiltonian with quadrupolc-quadrupolc correlation [87]. 

In Ref. [82], It was found that to a very large extent normal ordering (i.e., the 
sequence 0, 2, 4, etc.) of / in the yrast band can be kept when one changes the 
Hamiltonian from the realistic effective interaction to the TBRE Hamiltonian. This 
regularity is very interesting, but its origin has been discussed. 

Constraints on random interactions were studied by Kusnezov et al. [90] within 
the framework of the IBM and by Zhang et al. [92] within the geometric collective 
model. These works obtained some constraints on random interaction parameters 
of these models in order to produce global properties exhibited in the low-lying 
states of atomic nuclei. Within the framework of the shell model Horoi et al. found 
[94] that the B(E2) value obtained by random two-body interactions is too small 
in comparison to those obtained by the realistic effective interaction. It is difficult, 
however, to obtain the constraints on the TBRE in the shell model for heavy nuclei 
because of its huge dimensionality. 
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7 Summary 



The present subject was stimulated by the discovery of spin ground state (0 g.s.) 
dominance in the presence of random two-body interactions in 1998 [12]. Because 
this discovery is both surprising and interesting, it has sparked off a sudden interest 
in the origin of the g.s. dominance in the ground states of even-even nuclei. It 
also led to a number of other discoveries, for example, the generic vibrational and 
rotational band structure within the frameworks of the vibron model and the IBM. 

In Sec. 2 we established the notations and conventions. We reviewed very 
briefly the models of nuclear structure before we defined Hamiltonians for different 
systems involved in this paper. Then we presented an introduction to Monte Carlo 
samphngs, and defined the two-body random ensemble (TBRE). It was noted that 
the statistical patterns obtained by using other random two-body ensembles with an 
average being zero (such as the so-called random quasiparticle ensemble) are similar 
to those obtained by using the TBRE. 

In Sec. 3 we presented typical results of / g.s. probabilities, P{I), for various 
systems including fermions in a single-j shell, many-j shells, and d-, sp-, sd- and sdg- 
boson systems. In systems with an even number of fermions the g.s. probability is 
usually dominant with few exceptions, while for an odd number of fermions, where 
no / = state exists, the I — j {j is one of the angular momenta of single-particle 
states) g.s. probability is large. For systems with an even number of bosons the 

g.s. probability is large, while for an odd number of bosons the probability for 

1 — I {I is one of the spins of the bosons) g.s. is larger than for 7 = g.s. in many 
cases, although there may be 7 = states in these systems. The pattern of P(/)'s 
is sensitive to the displacement of the TBRE, except for fermions in a single-j shell 
and bosons with spin /. 

The parity distribution in the ground states of the TBRE Hamiltonian was found 
to be similar to that of realistic nuclei with mass number A larger than 120. The 
single-particle levels of these nuclei involve of both positive parity and negative 
parity. The observed parity of the ground states of these nuclei is always positive 
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for even-even nuclei, and otherwise it is positive or negative with about 50% for each. 
We showed that in the presence of the TBRE Hamiltonian, positive parity dominates 
in the ground states of even-even nuclei, despite the fact that the numbers of states 
with positive and negative parity are very close to each other. Because parity is a 
quantity which is much easier to handle than spin, we expect a sound understanding 
of parity distribution in the ground states obtained by using the TBRE Hamiltonian 
in the near future. 

In Sec. 3, we also showed that the odd-even staggering of binding energies arises 
from random two-body interactions. As for the effect of random interactions of 
higher rank, the discussion restricted to the sd bosons. According to Ref. [17], 
the features obtained from random Hamiltonians including three-body interactions 
do not change very much if the number of bosons is much larger than the rank of 
Hamiltonian. 

In Sec. 4 wc discussed some simple systems in which either the eigenvalues 
are linear in the two-body matrix elements or one can classify the two-body matrix 
elements in a simple way. Three techniques have been developed: the first technique 
is based on the geometry of eigenvalues [65] , and is applicable to systems in which the 
eigenvalues depend linearly on the two-body interactions. Such examples include 
d bosons and fermions in a j < 7/2 shell. The second technique is called the 
mean-field approach [51, 53], and is applicable to the cases where one can classify 
two-body interactions according to different geometric shapes of the systems. Such 
examples include sp bosons and sd bosons. The discussion of sp bosons based on 
random polynomials was also presented along this line. The third technique is called 
the empirical approach [56, 58]. Here one needs to know the lowest state and the 
highest state when one particular two-body matrix element is —1 and all others are 
zero. This approach can be applied to all simple systems discussed above. 

In Sec. 5 we discussed -P(/)'s of comphcated systems, for which one cannot 
classify the two-body matrix elements as done for the sp and sd bosons and the 
eigenvalues are not linear in the two-body interactions. In this case the empirical 
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approach was found to predict the -P(/)'s reasonably well. Here one needs the 
number of times for / to be ground state spin when one two-body matrix element is 
set to be —1 and others zero, and the same procedure is repeated for aU the two-body 
matrix elements. It was also demonstrated that the g.s. dominance in the presence 
of the TBRE Hamiltonian arises essentially from the two-body matrix elements 
which give 7 = ground states in this process. Differences between P(/max)'s of 
fermions and those of bosons were found and understood by using this empirical 
approach. 

For some systems such as four fermions in a single-j shell and four bosons with 
spin I, an argument was given that the g.s. dominance is partly due to the fact 
that there is only one non-zero eigenvalue for 7 = states when only one of the 
two-body matrix elements is switched on. 

In Sec. 5, we also reviewed an alluring but controversial argument of the 
g.s. dominance based on time reversal invariance of the Hamiltonian, the efforts to 
understand the g.s. dominance based on the large width of the distribution of 
eigenvalues for 7 = states, and an observation of a large pair-transfer fractional 
collectivity for g.s. of systems with particle numbers differing by two. 

Some by-products were obtained. For example, the ground state spin 7 of n 
fermions in a single-j shell was found to be equal to n when G2 — —1 and others 
zero {n is even). 

In Sec. 6, we first reviewed energy centroids of states with spin 7, denoted by 
Ej. The probability for Ej to be the lowest energy is large only when 7 ~ or 
7 ~ 7max- One thus can approximately divide the TBRE into two subsets, one which 
produces Ej^i^^^ as the lowest energy, and the other which produces Ei^i^^ as the 
lowest energy. The Ej averaged over each subset of the TBRE follows the 7(7 -|- 1) 
relation approximately. These features were explained in terms of the chaoticity of 
two-body coefficients of fractional parentage. The geometric chaoticity related to 
the regularities of energy centroids was also discussed by Zelevinsky and Volya in 
Ref. [9]. 
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In Sec. 6, we also showed that generic vibration and rotation arise in sp- and 
S(i-boson systems with TBRE Hamiltoniaiis. As for fcrmion systems, one needs 
additional requirements in order to contrive an ensemble that exhibits rotational 
behavior. To this end, two Hamiltonians have been suggested so far: One involves 
a displaced TBRE with an attractive average, the other involves the quadrupole- 
quadrupole correlation. 

In the third part of Sec. 6, we have seen that there is a large probabihty that 
the yrast states keep the "correct" ordering (i.e., 0, 2, 4 • • •) even for purely random 
interactions. In the fourth part, we discussed constraints on the TBRE Hamiltonian 
for the sd IBM [32] and geometric collective model [93] in order to reproduce the 
global properties exhibited by atomic nuclei. 

Parallel to the studies in the context of nuclear structure discussed in this paper, 
interesting features of randomly interacting quantum systems have been discussed, 
such as metallic clusters [18], quantum dots [95] and etc. [79, 96]. These works sug- 
gests that quantum systems or the dynamics of many-body systems in the presence 
of random interactions is an interdisciplinary and new field in theoretical physics. 

In conclusion, there are many interesting regularities for many-body systems in 
the presence of the two-body random ensemble, as we extensively discussed in this 
paper. Many of these regularities have been well explained theoretically. The 
g.s. dominance has been well confirmed by a large amount of numerical experi- 
ments. Many efforts have been devoted to understand why the g.s. dominance 
occurs. However, the underlying physical origin of the g.s. dominance is still to 
be explained in future. 

Acknowledgement We are also grateful to Drs. W. Bentz, G. Bertsch, R. 
Bijker, N. D. Dang, J. N. Ginocchio, V. K. B. Kota, and A. P. Zuker for their 
reading of this manuscript. We gratefully acknowledge interesting discussions with 
Drs. R. Bijker, R. F. Casten, S. Drozdz, A. Prank, J. N. Ginocchio, P. Van Isacker, 
C. W. Johnson, V. K. B. Kota, B. R. Mottelson, S. Pittel, O. Scholten, N. Shimizu, 
I. Talmi, V. Zelevinsky, and A. P. Zuker. The authors would hke to thank Dr. R. 



82 



Bijker for his allowing us to use Figures 18 and 19, which were taken from his papers. 



83 



TABLE I. g.s. probabilities for different random ensembles, as compared to 

the percentage of all states in the model spaces that have these quantum numbers. 
Data are taken from Table 1 of Ref. [13]. 



Nucleus 


TBRE 


RQE 


RQE-NP 


RQE-SPE 


1 = 












(total space) 




50% 


68% 


50% 


49% 


11.1 % 


220 


71% 


72% 


68% 


77% 


9.8% 


240 


55% 


66% 


51% 


78% 


11.1% 


^^Ca 


41% 


70% 


46% 


70% 


5.0% 


46Ca 


56% 


76% 


59% 


74% 


3.5% 


4«Ca 


58% 


72% 


53% 


71% 


2.9% 



TABLE IL Observed parity distribution in the ground states of atomic nuclei. 
We included all available data with mass number A larger than 120. The single- 
particle levels of the shell model space for these nuclei involve both positive and 
negative parity. The statistics here is based on the compilation in Ref. [60]. 



counts 


even-even 


odd-A 


odd-odd 


verified (+) 


361 


182 


68 


verified (-) 





164 


78 


tentative (+) 





109 


53 


tentative (-) 





109 


53 
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TABLE III. Probabilities for the ground states (in %) calculated by using the 
TBRE Hamiltonian to have positive parity. In brackets the number of neutrons and 



protons (iVp, N^) is given for each basis. 



basis A 












(0,4) (0,6) 
86.8% 86.2% 


(2,2) 
93.1% 


(2,4) 
81.8% 


(2,6) 
88.8% 






(2,3) (1,4) 
42.8% 38.6% 


(1,3) 
77.1% 


(0,5) 
45.0% 


(1,5) 
69.8% 


(6,1) 
38.4% 


(2,1) 
31.2% 


basis B 












(2,2) (2,4) 
72.7% 80.5% 


(4,2) 
81.0% 










(3,4) (3,3) 

'iZ.o/i / i.y/( 


(2,3) 
/ ZA/i 


(5.1) 

i-Z.y/i 


(3.2) 
■jy.i/i 


(4.1) 


(1.4) (5,0) 

Of""; {(''■/ II "1 ('/ 
lb A A 11.1/0 


basis C 












(2,2) (2,4) 
92.2% 81.1% 


(4,0) 
80.9% 


(6,0) 
82.4% 








(1,3) (1,5) 
73.0% 64.4% 


(2,3) 
52.0% 


(5,0) 
42.6% 


(4,1) 
56.5% 






basis D 












(2,2) (4,2) 
67.2% 76.1% 


(2,4) 
74.6% 


(0,6) 
83.0% 








(3,3) (3,2) 
54.5% 54.2% 


(2,3) 
54.0% 


(0,5) 
45.9% 
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TABLE IV. Probability for each state to be the ground state, and the distribution 
width of each eigenvalue for a j = 7/2 shell with four fermions. Each of the states is 
labeled uniquely by its angular momentum / and seniority number v. Probabilities 
of the row "TBRE" are obtained by 1000 runs of the TBRE Hamiltonian, and those 
of "predl." are obtained by calculating integrals such as Eq. (20) for O"*" state of 
n = 4, J = I case. The row "pred2." is obtained by using the empirical approach of 
Eq. (24) in Sec. 4.2. The row "exact" is obtained by using the geometry method 
proposed by Chau et al. [65]. The distribution width, gi(y)-i of each eigenvalue, is 
hsted in the last row. 



I{v) 


0(0) 


2(2) 


2(4) 


4(2) 


4(4) 


5(4) 


6(2) 


8(4) 


TBRE 


19.9% 


1.2% 


31.7% 


0.0% 


25.0% 


0.0% 


0.0% 


22.2% 


predl. 


18.19% 


0.89% 


33.25% 


0.00% 


22.96% 


0.00% 


0.02% 


24.15% 


pred2. 


14.3% 


0% 


28.6% 


0% 


28.6% 





0% 


28.6% 


exact 


18.33% 


1.06% 


33.22% 


0% 


23.17% 





0.05% 


24.16% 


9i(v) 


3.14 


3.25 


4.12 


3.45 


3.68 


3.62 


3.64 


4.22 



TABLE V. Angular momenta which give the largest (smallest) eigenvalues when 
Q = —1 and other q/ (/' ^ I) parameters arc for d boson systems. Here = 5, 
and A/"o =3, 0, 1, 2, 1, 0, for 6k,, 6k,+ 1, 6/t + 2, ■ ■ •, 6k + 5, respectively, k is a natural 
number. A/"/^^ = 2 for all n, and A/"2 = 5- Aq. This suggests periodical P(/)'s versus 
n for d bosons, according to Eq. (24) in Sec. 4.2. This table is obtained based on 
the reduction rule for C/(5) — > 0(3). 

n Co(min) C2(min) C2(max) C4(min) C4(max) 

6k 

6k+1 2 2 

6k+2 2 

6k+3 2 

6k+4: 2 

6/V+5 2 2 



-'max -'max 

T T 2 

max max ^ 

T T 2 

max max 

-^max -^max 

T T 2 

max ^ max ^ 

I I 

-'^ max -'^iiiax — 
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TABLE VI. Angular momenta which give the lowest eigenvalues when Gj = — 1 

and all other two-body matrix elements arc zero for four fermions in a single-j shell. 



2') 








(la 


do 

^8 




G\2 


Gi4 


G16 


G18 




G22 


G24 


c 


c 

^28 




7 


n 
u 


A 


9 


Q 




























n 

y 


U 


4 


U 


U 


1 
























11 





4 





4 


8 


16 






















13 





4 





2 


2 


12 


20 




















15 





4 





2 








16 


24 


















17 





4 


6 





4 


2 





20 


28 
















19 





4 


8 





2 


8 


2 


16 


24 


32 














21 





4 


8 





2 











20 


28 


36 












23 





4 


8 





2 





10 


2 





24 


32 


40 










25 





4 


8 





2 


4 


8 


10 


6 





28 


36 


44 








27 





4 


8 





2 


4 


2 








4 


20 


32 


40 


48 






29 





4 


8 








2 


6 


8 


12 


8 





24 


36 


44 


52 




31 





4 


8 








2 





8 


14 


16 


6 





32 


40 


48 


56 



TABLE VIL Same as Table V for sd-boson systems. The angular momentum 
I corresponding to Csdsd = — 1 is omitted because it always presents degenerate 
levels for states of many 7's. The one-body parameter is omitted because we are 
interested in P(/)'s in the presence of the TBRE Hamiltonian. 



n 


^ssss 




(^ssdd 


Co 


C2 


C4 


6 

















-^max 


7 











2 


2 


-^max 


8 














2 


-^max 


9 











2 





-^max 


10 














2 


-^max 


11 











2 


2 


-^max 


12 

















-^max 


13 











2 


2 


-^max 


14 














2 


-^max 


15 











2 





-^max 


16 














2 


-^max 
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TABLE VIII. Coefficients a/ and V{I) for four fermions in a j = | shell. Bold 

font is used for the largest aj, and italic for the smallest aj for a given J. Probabil- 
ities in the column "predl." are obtained by integrals similar to Eq. (7) in Ref. [47], 
and those in the column "pred2." are obtained by the empirical formula given in Eq. 
(24). The P(J)'s in the last column "TBRE" (in %) are obtained by diagonalizing 
the TBRE Hamiltonian for 1000 runs. We take both the smallest and the largest 
aj when counting J\fj. 



I 


Go 


G2 


G4 


Ge 


Gs 


predl. (%) 


pred2.(%) 


TBRE 





0.80 


0.35 


1.74 


2.11 


1.01 


11.97 


11.1 


10.2 


2 


0.30 


1.39 


1.45 


1.29 


1.56 


14.51 


22.2 


15.4 


3 


0.00 


0.36 


2.28 


2.63 


0.71 


28.17 


33.3 


28.9 


4 


0.20 


1.07 


1.38 


1.91 


1.44 


1.74 





1.7 


5 


0.00 


1.00 


1.59 


1.84 


1.57 


0.30 


06 


0.6 


6 


0.20 


0.79 


1.50 


1.58 


1.93 


0.22 





0.3 


7 


0.00 


1.20 


1.09 


1.40 


2.31 


3.44 





3.2 


8 


0.30 


0.48 


1.05 


1.82 


2.36 


0.03 








9 


0.00 


0.17 


1.33 


2.12 


2.38 


0.01 








10 


0.00 


0.70 


0.69 


1.41 


3.21 


6.76 





8.7 


12 


0.00 


0.00 


0.52 


1.69 


3.78 


32.64 


33.3 


31.0 
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Figure captions: 

FIG. 1 Probabilities of /+ ground states for different j shells with four fermions. 
All probabilities are obtained from 1000 runs of the TBRE Hamiltonian. One sees 
that P{0) periodically staggers with the value of j at an interval Sj = 3. 

FIG. 2 Same as Fig. 1 for n=5. 

FIG. 3 P(0)'s of n = 4, 6 and -P(j)'s of n = 5, 7 fermions in a single-j shell. 
They stagger synchronously at 5j — 3 when j is small and seem to saturate when j 
becomes large. 

FIG. 4 I g.s. probabihties vs. I of four bosons with spin I. The results are 
obtained by 1000 runs of the TBRE Hamiltonian. The P(0) staggers with spin I at 
an interval Si — 3 (similar to the P(0) of four fermions in a single-j shell, refer to 
Fig. 1), and that the P(/max) is very large. 

FIG. 5 7 = and I — I g.s. probabilities versus n ior I — A and 6. One sees 
that P{0) is usually smaller than the corresponding P{1) when n is an odd number, 
indicating that the g.s. dominance might be associated with an odd-even effect of 
boson number. 

FIG. 6 I g.s. probabilities for d bosons, with boson number n ranging from 
4 to 44. Only states with I —0, 2, and /max = 2n are possible as ground states. 
g.s., 2 g.s. and /max = 2n g.s. probabilities are periodically close to 0, 20%, 40% or 
60%. P(0) ~ when = 6k ± 1. The predicted P(/)'s (open squares) are well 
consistent with those (sohd squares) obtained by using the TBRE Hamiltonian. 
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FIG. 7 Polygons corresponding to a system of five d bosons. Each state is 
represented by a dot. The dots inside tlie polygon never come to the ground, and 
the g.s. probability of each state of the vertex is determined by Eq. (31). Here 6i 
(corresponding to the I{v) = 2(1) state) and 9^, (corresponding to the I{v) = 2(5) 
state) lead P(2) ~ 60%, 64 leads to P(/max) ~ 40%, and 62 (corresponding to the 
I{v) = 0(3) state) leads to P(0) ~ 4%. 

FIG. 8 -P(O), -P(l) and P(n) for sp bosons. The results are obtained by 1000 
runs of the TBRE Hamiltonian. 

FIG. 9 P(0)'s of fermions in a single-j shell. Sohd squares are obtained by 

1000 runs of the TBRE Hamiltonian. The open squares are predicted P(0)'s. a) 
n = 4; b) n = 6. Solid triangles are obtained from the empirical formula of Eq. 
(36). 

FIG. 10 Comparison of P(/)'s obtained by diagonalizing the TBRE Hamilto- 
nian with those predicted by the empirical rule of Eq. (35). Here we show fermions 
in two-j shells with (71,72) = (|i §)• =4, 5, 6, 7 in a), b), c) and d), respectively. 
Solid squares are obtained by 1000 runs of the TBRE Hamiltonian and open squares 
are the predicted values. 

FIG. 11 The P(0), P(2) and P(/max) of sii-boson systems. Solid symbols 
are P(/)'s obtained from 1000 runs of the TBRE Hamiltonian. Open symbols are 
P(/)'s predicted by the empirical approach introduced in Eq. (35). Only / = 0, 2, 
-^max g-S. probabilities are included. All other P(/)'s obtained by diagonalizing the 
TBRE Hamiltonian are close to zero, and the predicted P(/)'s are also zero. 
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FIG. 12 a) The /max g-s. probability obtaining by fixing G'i6=±l (Jmax=16) 
and all other Gj being a TBRE multiplied by e; b) g.s. probabilities obtained 
by fixing G'o=±l and all other Gj being the TBRE Hamiltonian multiplied by e. In 
this figure, j = y n = 4. 

FIG. 13 Regularities of /maxg-S- probabilities. The squares are predicted by 
the 1/A^ relation whereas all other results were obtained by diagonalizing the TBRE 
Hamiltonian. (a) Fermions in a single-j shell; (b) /j^j^^^ g.s. probabilities for fermions 
in two-j shells; (c) bosons with spin /. The I'^^^ is defined by ImaxUi) and /max(i2)- 
One sees that the 1/A^ relation applies well to the fermions in a single-j shell but 
underestimates the P{Imax} for bosons with spin I when I becomes large. The 1/N 
relation predicts reasonably the lower limit of -P(/max) fermions in two-j shells. 
Refer to the text for details. 

FIG. 14 Seniority distribution in the angular momentum / = ground states. 
No bias of low seniority is observed in these systems with four and six fermions in 
a single-j shell. 

FIG. 15 Typical results for 7^(1) 's: a) single-j (j = y) with four fermions, 

b) two-j shell (2ji,2j2) = (11,7) with six fermions, c) 6-sdg bosons, d) single-j 
U — 2) with five fermions. 

FIG. 16 {Ei)^-aia vs. /(/ + 1). (i?/)inin's are obtained by averaging over Ej 
with the requirement that Eir^j^^^.^^ be the lowest energy. The results are obtained by 
applying 1000 sets of the TBRE Hamiltonian. a), twenty d bosons; b) ten sd bosons; 

c) single-j shell with j=17 and n — A, and d) a two-j shell with ji = 5/2, j 2 — 7/2, 
and n —A. The quantity {Ei)ja.m has the same behavior except for a difference of 
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sign. Refer to the text for details. 



FIG. 17 Correlation between ^ and j = ySil- ^ - l-42yEjF for sd 
bosons, and fermions in a single- j shell. This correlation is shifted very shghtly to 
the right for fermions in many-j shells and sdg bosons. 

FIG. 18 Distribution of R of sd bosons in the presence of random two-body 
interactions with n =3 (dash-dotted), 6 (dotted), lO(dashed) and 16 (solid). 

FIG. 19 Correlation between R and 4. — p4r( for sixteen sd bosons with one- 

body and two-body random interactions. It is seen that (i?, p/^o'o+^n+! ) concentrate 
at two points: (2.0, 2.0) which is characteristic for vibrational motion and (3.3, y) 
which is characteristic for rotational motion in the large n limit of the IBM. 

FIG. 20 The distribution of it! values for six identical nucleons interacting by 
the TBRE Hamiltonian in a) the sd shell, b) the p/ shell, and c) the sdg shell. 
The calculations were done within the SD-^qat truncated subspace. 

FIG. 21 The distribution of R values for six identical nucleons in the sd 
shell interacting by a random Hamiltonian containing monopole pairing, quadrupole 
pairing, and quadrupole-quadrupole forces. The strength k of the quadrupole- 
quadrupole interaction is multiplied by a factor e to assess the importance of the 
quadrupole-quadrupole interaction. The calculations were done within the 5'D-pair 
truncated subspace. 

FIG. 22 The distribution of it! values for six identical nucleons interacting 
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by a random Hamiltonian containing monopole pairing, quadrupolc pairing, and 
quadrupolc-quadrupolc forces for a) the pf shell, b) the sdg shell, c) the pfh shell, 
and d) the sdgi shell. The calculations were done within the SD-pa-ir truncated 
subspace, and e is fixed to be 1.0. 

FIG. 23 Correlation between ratios of E2 transition rates and R for the same 
calculations as in Fig. 22d . The inserts in panels a), b) and d) focus on critical 
regions. 
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Appendix 



List of mathematical notations used in this paper. 



A mass number of an atomic nucleus 

Z proton number of an atomic nucleus 

N neutron number of an atomic nucleus 

(Nn) valuence proton (neutron) number outside a closed shell 

n number of valence particles 

j angular momentum of a single — particle state 

I intrinsic spin for a boson or the orbital component of j 

m z — component of j 

t isospin for a nucleon 

rut z — component of t 

T total isospin of two or more nucleons 

Mt z — component of T 

J total angular momentum of two fermions or two bosons 

Mj z — component of J 

p{x) distribution function of x 

ejmt single — particle energy of valence protons or neutrons [rrit — ^-) 

(cp) single d (p) boson energy 
GjT{jij2-ijzji) two — body matrix elements of fermions 

Gj abbreviation of GjT{jij2-ijz3i) for single— j fermions 

Gl abbreviation of GLrihh, hh) for spin — I bosons 

TBRE abbreviation for "two — body random ensemble" 

N number of independent two — body matrix elements 

Nm 27V - 1 
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H[ number of times that spin / states appear either as the ground 

state or the highest state when one of Gjt(jiJ2J3J4)'s is — 1 
and others are zero 

Hi number of times that spin I states appear as the ground state for 

one of G jT{iii23?,34)'^ to be — 1 and others to be zero 
K parameter for quadrupole — quadrupole interaction 

/ total angular momentum for a state of many — body systems 

-^max maximum of / 

-^min minimum of I 

V seniority number for fermions in a single— j shell 

or number of d bosons not paired to spin zero 
(5 additional quantum numbers (except /) to specify a state 

S (D) fermion pairs with spin zero (two) 

P{I) probability that the ground state has angular momentum / 

Ej energy centroid of spin I states 

V{I) probabihty that Ej is the lowest among all Eji 

Dj number of the angular momentum / states 

Dj'^ Dj for fermions in a single— j shell 

Df^ Dj for bosons with spin I 

Eq^''^ the non — zero eigenvalue for / = of four bosons with spin / 

and Gl = Sll' 

Eq^^ the non — zero eigenvalue for 7 = of four fermions in 

a single— j shell and Gj — —5jj> 
ai width of spin / states, defined by a] = {{H — Ej)"^) / Dj 

9i{v) width of the /(v) state, defined by 
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a 



( ) min 
{El) max 
J 

I)) 

R 

/transfer 



^ 7 
X/ i'^Hv)) ' fermions in j < - shell 



width of spin I states , defined by gj — {H^)/Di 



ensemble average of Ei with £'/^/^.^being the lowest 



ensemble average of Et with -E"/^/ being the lowest 

>^ (-£'7) min 



coefficient obtained by (-E7)inin = "^^"[^^ 



or by {Ei)m^ ^ Eo - 



2J 



normalized state of | ) 

excited energy of the first state of an even — even nucleus 
E4+/E4+ 

pair — transfer fractional collectivity 
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